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1 Introduction 



In this paper we review the formahsm of free fermions introduced by the Kyoto school 
and used for construction of tau-functions of classical integrable hierarchies of 
nonlinear PDE's such as Kadomtsev-Petviashvili (KP), modified Kadomtsev-Petviashvili 
(MKP), 2D Toda lattice (2DTL) and their multicomponent versions. The operator ap- 
proach to the integrable hierarchies appears to be an extremely useful tool in dealing 
with classical integrability because the bilinear equations of the hierarchies acquire the 
transparent meaning as the standard relations for correlation functions of free fermions 
like the Wick theorem. The r-function itself is a vacuum expectation value of special 
operators from the infinite-dimensional Clifford algebra generated by free fermions ipn, 

We give a detailed and self-contained derivaton of the key properties and relations: 
group-like properties of the normally ordered exponents, transformations between dif- 
ferent normal orderings, the bilinear relations, the generalized Wick theorem and the 
boson-fermion correspondence, paying attention to particular cases and examples. Our 
motivation comes from the feeling that their precise formulations and exhaustive proofs 
seem to be missing in readily available literature. In view of this, we believe that our 
review might be of interest for those who are going to use the operator methods in their 
own research work. 

Another motivation is to collect together various examples of r-functions that play an 
important role in the theory and/or applications. Among them are familiar ones (char- 
acters, quasi-polynomial and multi-soliton r-functions, partition functions of different 
matrix models) as well as less known examples which were addressed in the literature 
only recently (such as expectation values with non-standard time evolution and Hur- 
witz partition functions). In each case we present explicit fermionic realizations of the 
r-functions. 

In order not to overload the paper, we restrict ourselves by one-component charged 
fermions and leave aside integrable hierarchies of the types other than Aq^- Some other 
related interesting topics, such as fermionic representation of the generalized Kontsevich 
model and the melting crystal model, remain out of the scope of this review because 
of the space limitations. Various aspects of the free fermionic approach to integrable 
systems are also discussed in the books [5], chapter 9], P, chapter 14] and reviews [3 
section 4], [HI section 1]. 

The organization of the paper is clear from the table of contents. It should be em- 
phasized that central in our presentation is the basic bilinear condition (BBC) 

i^kG ®rkG = Y. Gi^u ® Grk (i.i) 

k k 

which defines an important class of elements G of the Clifford algebra. Exponents (or 
normally ordered exponents) of bilinear expressions in the fermionic operators ipn-, V'n ^ire 
well-known solutions to the BBC. There are also solutions which can not be represented 
in this form - for example, any linear combinations of the tpn^ (or V'n's)- Correlation 
functions of fermions with insertions of elements G that obey the BBC have some special 
properties. In particular, the generalized Wick theorem holds for them. Taking matrix el- 
ements of the operator relation f 1 1.1 1) between some properly chosen states in the fermionic 
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Fock space, one can derive the bilinear identity generating various bilinear equations of 
the Hirota type for the r-functions. The general r-function for one-component fermions 
is usually defined as the expectation value of the form 



where t± = {t±k}'^^i and n = to E 1j are time variables, J± are current operators and G is 
any solution to the BBC with zero charge. An extension of this definition to solutions of 
the BBC with non-zero charge is straightforward. By freezing some of the time variables 
at particular values (say, t_ = Oort_ — n — 0) one obtains r-functions of the MKP and 
KP hierarchies respectively. 
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2 Free fermions 



2.1 The algebra of fermionic operators 

Let ipnj'ipn^ n e Z, be free fermionic operators with usual anticommutation relations 

[ijjn, '0m]+ = [C -0^]+ = 0> bl^n, "0^]+ = <^mn- (2.1) 

They generate an infinite dimensional Clifford algebra. We also use their generating 
series 

From the fact that the anticommutator of any linear combinations of the fermionic 
operators is a number it follows that the commutator of any bilinear expressions in -0^ 
and is again bilinear in ipn and For example, 

We see that the commutation law for the expressions iprni^n same as for the matrices 
Emn with matrix elements {E^n)ij — ^im^jn which are generators of the algebra gl{oo) of 
infinite-size matrices with only finite (but arbitrary) number of non-zero elements. More 
generally, consider the bilinear expression 

XA = Y,Aj^ljiij* (2.3) 
with some matrix A, then [Xa, Xb] — Xia,b] and 

i i 

In order to find the adjoint action of the operator Xa on fermions, we use the well 
known formula 

e^Be-^ = B + [A,B] + ^[A, [A, B]] + . . . (2.4) 
valid for any two operators A, B. We get 

e^^V'ne-^^ = E ' e^^Ce-^^ = J2{R-\iij*, (2.5) 

i i 

where R^^ — e"^^ = 1 ±^4-1- ^A^ + We see that exponents of the expressions bilinear 

in ipn and possess a rather specific property: the result of their adjoint action on the 
form linear in fermions is again linear. 

2.2 Dirac vacua and excited states 

Next, we introduce a vacuum state |0) which is a "Dirac sea" where all negative mode 
states are empty and all positive ones are occupied: 

i^n\0) = 0, n<0; dO) = 0, n > 0. 
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(For brevity, we call indices n > positive.) Similarly, the dual vacuum state has the 
properties 

(0|V': = 0, n<0; (0|V'n = 0, n > 0. 

With respect to the vacuum |0), the operators ipn with n < and ip* with n > are 
annihilation operators while the operators ip* with n < and ipn with n > are creation 
operators. 

We also need "shifted" Dirac vacua \n) and {n\ defined as 

V'n-i ■ ■ ■ V'lV'o |0) , n>0 
\n) = { (2.6) 
C---V'-2V'-i|0), n<0 

{n\ = { (2.7) 

(0| llj-llp-2 ■■■ipn, n<0 

In particular, we have 

ipm 1^) =0, m < n; ipm 1^) =0; m > n, (2.8) 

(n| iprri — 0, m> n; {n\ -0^ = 0, m < n. (2.9) 

It is also clear that 

iPn\n) = \n + l) , rn\n+l) = \n), (2.10) 

(n+l|Vn = (ri| = + (2.11) 

Excited states (over the vacuum |0)) are obtained by filling some empty states (acting 
by the operators ip*) and creating some holes (acting by the ip/s). A particle carries 
the charge —1 while a hole carries the charge +1, so any state with a definite number 
of particles and holes has the definite charge. Let us introduce a convenient basis of 
states with definite charge in the fermionic Fock space Hp- The basis states \X,n) are 
parametrized by n and Young diagrams A in the following way. Given a Young diagram 
A = (Ai, . . . , Xi) with i = i{X) nonzero rows, let (a|/3) = (cti, . . . , aci(\)\/3i, . . . , (3d(\)) be 
the Frobenius notation for the diagram A (see Appendix A). Here d{X) is the number 
of boxes in the main diagonal and ai = Xi — i, /3i = X[ — i, where A' is the transposed 
(reflected about the main diagonal) diagram A. Then 

|A, n) := V;-/3i-l • • • V'n-/3,(;,)-l V'n+a,(,) • • • i^n+ai \n) , 

(2.12) 

(A, n\ := {n\ ^/j*^^^ . . . i'n-^aw-^ ■ ■ ■ V'n-/3i-l- 

The state \X,n) has the charge n with respect to the vacuum state |0). For the empty 
diagram (0,n| = {n\, |0,n) = \n). 

The states |A, n) can be constructed from a vacuum in another, equivalent way which 
is sometimes more convenient. 
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Proposition 2.1 We have: 

\X,n) = (-l)^(^) ^pn+\i~li'n+\2-2 ■ ■ ■ i'n+Xe-eln - i) 



and 



(A, n\ = (-l)^W {n - i\ rn+x,-e ■ ■ ■ C+a.-2^:+a,-i 
I A, n) = C-A^+i • • • C-a;„+™-i + ^) 



(A,n| = (-1)1^1 ^^^^ {n + m\'l/jri-\'^+m-l---i'n-X'^+li'n-\[ 

where i = i{X) = X'l, m = Xi and 

d(A) 

6(A) = 5^ (A + 1). 



(2.13) 



(2.14) 



(2.15) 



i=l 



This statement is basically a rephrasing of Proposition (1.7) from Macdonald's book [9], 
with an additional input of commutaton relations for fermionic operators. The proof is 
a simple combinatorial exercise. We give an idea how to prove ( I2.13p . Denote the r.h.s. 
by \X,n) , then we can write 

|A, n)' = V^„+Ai-l . . . ^n+Xa-d ^n+A,+i-(d+l) • • • V^n+A,-£ ^n-i ■ ■ ■ 1^) 

e-d {e-d)+d 

where d = d{X) is the size of the maximal square of the diagram. Since Aj — z = at 
i = 1, . . . ,d, we see that the first d ■^/'-operators are exactly those standing in f l2.12p . The 
£ — d ■j/'-operators from the next group have indices which are strictly less than n, and each 
^/'-operator ipn-i from this group has its ?/'*-partner ip^-i among the '?/'*-operators from the 
third group. These pairs ipn-ii^n-i cancel when acting to the vacuum \n) and we are left 
with just d '?/'*-operators which are exactly those standing in fl2.12p . (It may be helpful 
to note that the Frobenius notation corresponds to Aj's as follows: {Aj — i}^!^'' = {aijfl'^'*, 

- ^^}^dlx)+l = {1, 2, . . . , /3i + 1} \ {A + 1}£M Therefore, |A, n)' differs from |A, n) 
only by a sign. 

Remark. Note that if one formally adds some "zero parts" with Aj = to the partition 
A, then the state fl2.13p remains unchanged. 

From the definition of the vacuum states \n) it is obvious that for any such state 
and arbitrary /c G Z either ipk or ipl is an annihilation operator: either ipk\n) = or 
t/'^ |n) = 0. From definition fl2.12p it is also obvious that the same is true for any basis 
state: 



and 



where /„ a is the set 



'n,A 



iJk \X,n) = 0, 




V^fc|A,n) = 0, 


k i In.X 


ipk \X,n) ^ 0, 


^ ^ in,X\ 




k e In,X 


(A,n| V'fc = 0, 


k ^ -^71, A) 


(A,n|7/>* = 0, 


k e In,X 


(A,n|?/'fc ^ 0, 


k e In,X] 


(A,n|V'^^0, 


k i In.X 



(2.16) 



(2.17) 



A; G Z 



<n,fc^{n-A-l}2)}u{n + a.};ri 



,d(A) 
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2.3 The vacuum expectation values 



The vacuum expectation value (0| . . . |0) is a Hermitian linear form on the Clifford algebra 
fixed by 

(0|0) = 1. 

Then, from the commutation relations (12. ip and definitions of the "shifted" Dirac vacua 
fl2.6p . (12. 7p it follows that (n|n) = 1 for any n. Bilinear combinations of fermions satisfy 
the properties (n| ipiipj \n) = {n\ \n) = for all i,j and 

{n\ ipiip* \n) = 6ij for j < n, {n\ ipiip* \n) = for j > n. 

The expectation value of any operator with non-zero charge is zero. The basis vectors 
f l2.12p are orthonormal with respect to the scalar product induced by the expectation 
value: 

(A,n| /i,m) = 5mn5\i,. 

This can be directly seen by taking the scalar product of the vectors of the type (12.121) 
and moving the operators ipn-Pi-i to the right, taking into account that the sequences 
«!, a2, . . . , ctd and /3i, /32, ■ ■ ■ , /3d are strictly increasing. 

From orthonormality of the basis states it follows that if {X,n \U) = (respectively 
{U\ X,n) = 0) for all A and n then \U) = (respectively {U\ = 0). Moreover, if a Clifford 
algebra element X kills all basis states, that is X\ X,n) = (or {X,n\X = 0) for all A 
and n, then X = 0. 



Proposition 2.2 The coefficients {X,n \U) (respectively {U\ X,n)) uniquely determine 
the state \U) (respectively {U\): 

|f/)=5^|A,n) {X,n\U) 

Similarly, the coefficients {X,n \X\ fi,m) uniquely determine the Clifford algebra element 
X: 

X= \X,n) {X,n\ X \fi,m) {fi,m\ . 

n,\;m,fi 

The explicit form of the operators |A,?t,) {fi,m\ in terms of fermions is given below in 
section 12. 8[ 

In general, expectation values of products of fermionic operators are given by the 
Wick theorem. Let Vi = Vijipj be linear combinations of ipf^ only and w* = w*jipj 
be linear combinations of ip^s only, then the standard Wick theorem states that 

(n| f 1 . . . VmW^ . . .wl\n) = det (n| ViW* \n) , 

i,j=l,...,m ■' 

(n| wl . . . w^Vm ■ ■ - Viln) = det (n| w*Vj \n) . 

i,j=l,...,m 

The Wick theorem can be proved by induction. We will not give the proof here because 
the proof of a more general statement will be presented below. 
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For the generating series ip{z), ip*{C) we have: 

I ncmz) \n) = J2 c'z' (n| r.i'k \n) = Y.^z/cf = 

j,k k>n 



in 



(assuming that \(\ > \z\) and 

I mriO \n) = J2 C'z' (n| \n) = Y.{z/Cf = 

j,k k<n 



in 



(assuming that \z\ > \(\). More generally, 

c 



(n 



riCi) ■ ..riCmmzm) . . .^(^i) i^^) = Hizi/cr ■ det _ 



Uizi-Zi') n (0 -Cj' 

i<i' j>j' 



(2.18) 



(assuming that > \zj\) and 

(n| 7A(^i) . . . ^PizmWiCm) . . . riCi) \n) = Hizi/CiT ■ det 

1=1 "'^ 

Uiz.-z,) n(o-O') 

(where \zi\ > \Cj\)- One can also calculate 

m m—l 
i<i' j>j' 



(2.19) 



m—l 



(n+/| riCm-i) . . . riCiMzi) . . . ^{zm) \n) = '-^^—^-^ n zi n c,V" 

r=l s=l 

(2.20) 

(where > \zj\) as a limiting case of f l2.18p by tending subsequently Cm, Cm-i • • • , Cm-i+i 
to infinity. In particular, 

m 

{n + m| 7/'(zi) . . . ^{zm) \n) = l[{zi - zj) J] 4. (2.21) 

i<j k=l 

In a similar way, we get from ( 12.19P 

m m—l 



' Zj Zji 



m m—l 



n(c -oiif 

{n-l\ ^(z„_0 . . . ^(zi)^*(Ci) . . . ^*(Cm) \n) = '^^^^;^T^ H ^^"11 ^fc"' 

n n(-^ - c^) 

r=l s=l 

(2.22) 
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(where \zi\ > and 

m 

(n-m| ^*(Ci) . . . riCm) \n) = 1[{Q - Q) J] C^^ (2.23) 

i<j k=l 

2.4 Normal ordering 

Once the vacuum state is fixed, a useful notion is the normal ordering of operators. The 
normal ordering ;(...); with respect to the Dirac vacuum |0) is defined as follows: all 
annihilation operators are moved to the right and all creation operators are moved to the 
left taking into account that the factor (—1) appears each time two neighboring fermionic 
operators exchange their positions. For example: li'li'i', = —ipi'ipi, I'^-iV'oI = ~'4'o4'-i, 
',ip2'ipii'ii'l2'. = '4^2'ipii'-2i'h 6tc. We also note the identity 

The normally ordered expressions are always well-defined when acting to the correspond- 
ing basis states of the fermionic Fock space Tip. 

Under the sign of normal ordering, all fermionic operators ipj and ip* anticommute. 
In other words, it is wrong to use the commutation relations of the Clifford algebra under 
the sign of normal ordering, i.e., for example, [ipi'^i'. 7^ l(l ~ V'l'^i)!- 

Using the normal ordering, one can introduce the charge operator Q as 

This operator counts the charge of the state: Q\X,n) = n\X,n) and so {fi,m\Q \ X,n) = 
nSnmSxn (note that without normal ordering this matrix element would be ill-defined!). 
The charge operator has the commutation relations [Q,ipn] = i^n, [Qyi^n] = ~4'n which 
mean that ipn, V'n have charges ±1. More generally, we say that a Clifford algebra element 
X has definite charge q if [Q,X] = qX. 

In a similar way, one may define the normal ordering with respect to any vacuum state. 
The general rule is that the annihilation operators are moved to the right and creation 
operators to the left (with the appropriate sign factor). With respect to the vacuum \n), 
the annihilation operators are ip* with j > n and ipj with j < n while ip* with j < n and 

with j > n are creation operators. We denote the corresponding normal ordering by 
',{■ ■ ■)',n (^^ ^^^^ notation ;(...); = ',{■ ■ ■)',q)- One can go even further and define the 
normal ordering with respect to any basis state 02.121) : the operators i/jk with k G In,\ and 
ipl with k ^ In,x are annihilation operators for the state |A,n), while all other fermions 
are creation operators, see (12.161) and (12.171) . 

One may also consider the bare vacuum \oo) which is the absolutely empty state. 
With respect to the bare vacuum, all ijj/s are annihilation operators while all ip^s are 
creation operators. Let us denote the corresponding normal ordering by . (for- 
mally, ^ (...) = l{. . Oloo)- ■'^^ ^^^^ normal ordering, all iIj*^s are moved to the left and 
all V^'s to the right, with taking into account the sign factor appearing each time one 
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operator is permuted with another. For example: ^^j^"^" x ~ '^m'^n, x^n'^/'mx ~ ~'^m'^n 
and 

I exp(^^ BikiJ-tpk^ X = 1 + Bik'ipii^k + ^ X] BikBek''ip*ip*,'ipk'i^k + ■■■ (2.25) 

The definition of the normal ordering is closely related to the definition of the vac- 
uum expectation value: ',4'k4'i'.n ~ '^I'^i ~ {'^I'^l'^i l"^)- More generally, for any linear 
combinations /o, /i, . . . , /m of the fermion operators ipi, ip* we have the recursive formula 

m 

/o:/i/2 • • • /m:„ = :/o/i/2 • • • + E^-^)'"' M- :/i/2 . . . | . . . (2.26) 

where ^- means that this factor should be omitted. This recursive relation allows one 
to express the normally ordered monomials with arbitrary number of fermions as linear 
combinations of monomials without normal ordering and vice versa. For example: 

'.foflf2'.n = /0/1/2 — (/l/2)„ /o + (/o/2)„ fl — (/o/l)„ /2, 
',foflf2f-i',n ^ /0/1/2/3 — {f2f-i)n /o/l + (/l/s)^ /0/2 " (/l/2)„ /o/s 
~ (/o/l)„ 72/3 + (/o/2)„ /1/3 — {fof-i)n /1/2 
+ (/o/l)n {hf3)n ~ (/o/2)„ (/l/3)„ + (/o/3)„ (/l/2)„ , 

where := (n| /Jj |n). 

2.5 Group-like elements 
2.5.1 Group elements 

Bilinear combinations bmni'mi^n of the fermions, with certain conditions on the 
matrix 6 = {bmn), generate an infinite-dimensional Lie algebra [2]. Exponentiating these 
expressions, one obtains an infinite dimensional group (a version of GL{oo)). Elements 
of this group can be represented in the form 

G' = exp(5^ feifcVr^fc). (2.27) 

i,kez 

The inverse element can be written in the same way with the matrix {—bmn)- 

As it was already stated, the group elements of the form (12.271) obey a rather special 
property that the adjoint action of such elements preserves the linear space spanned by 
the fermion operators ipn and the space spanned by More precisely, we have: 

I I 

or 

= Yl ^'"^^ ^ ' i^nG = J2 RnlG^l , (2.28) 

I I 
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where the matrix R = (Rni) of the induced hnear transformation is given hjR = e^. The 
product of two group elements is also a group element: 

exp ( Kkrii^k) exp ( kki':iJk) = exp ( 5^ (2.29) 

j.fceZ i,k£Z i,k£Z 

with e^'e^ = e^" . Clearly, multiplication of G of the form fl2.27p by any constant number 
preserves the characteristic property fl2.28p . 

We refer to the transformation fl2.28p as rotation of the fermionic operators with the 
rotation matrix R. From the fact that the center of the Clifford algebra is formed by 
numbers c [10] it follows that two group elements G, G' with the same rotation matrix 
R may differ by a scalar factor only: G' = cG. In particular, sometimes it is more 
convenient to consider the group elements with exponentiated normally ordered bilinear 
fermionic combinations, for example: 

exp ( X] feifc : i^ii'k : ) = Cfc exp ^ biki^*iJk) , 

where Cb = exp (- X]fc>o hk) ■ 

Note also that the charge operator Q f l2.24p commutes with any group element. 

2.5.2 Normally ordered exponents 

Let us note, following the works of the Kyoto school [21 [10], that the group elements 
can be equivalently represented as normally ordered exponents of bilinear forms. For 
example, one can directly prove that the element G = ^ ^Bikip'ipk x (\^qyq below 
summation over repeated indices is implied) satisfies the first commutation relation fl2.28p 

with Rin = 6in + Bin. 

= ^* XeB,.^*V. X + Ba,nra,+Ba,nBa,^A:,ra,i^b, + ^Ba,nBa,^,Ba,b,^^^^^ + • • • 

— ,/;* Xp-BifcVrV'fc X I D ^1,* X pBiktp*tpk X — (J) _|_ R X pBiktp*il)k X 

— Ynx^ x'-^an^^ax^ x ~ K'^an ^ J->an ) Y a x^ x' 

The second commutation relation fl2.28p can be proved in the same way. Moreover, any 
solution to ( I2.28P can be represented by a normally ordered exponent ^ e^iki't'^k x _ j^^ f^^j^ 
be also checked, in a similar way, that 

exp (foifcV^r^fc) = I exp ((e'' - I'U'^^l^k) I , (2.30) 

where I is the unity matrix. The composition law is given by 

I exp(5:,V^*^,) I I exp(5,,^;^,) ^ = ^ exp ((5 + fi'+M),,V^;^,) ^ (2.31) 

which directly follows from the composition law f l2.29p and the formula B = ~ I. 

Let us prove another useful formula which allows one to represent a group element as 
a normally ordered exponent with respect to different vacua: 
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Proposition 2.3 Whenever both sides are well-defined the following is true 

- exp(^B,kiJ*i^k) ^ = det (J+P+5) ; exp (^i^Vr^fc) I , (2-32) 

or, equivalently, 

: exp(^AikiJ*^k) : = det{I-P+A) - exp(^Bik^*^k) • (2.33) 

Here P+ is the projector on the positive mode space ({P^)ik = Sik for i,k > and 
otherwise) and the matrices A,B are connected by 

B-A = AP+B, I.e., B = {I-AP+)-^A or A = B{I+P+B)-\ (2.34) 



The idea of the proof is as in [TU] . First we notice that we can write the normally ordered 
expression •e^^''^^'^'' I as the composition of three operators: 

Gi G2 G3 

where the repeated indices a, 6 (a, 6) in the r.h.s. imply summation over non-negative 
(respectively, negative) integers (summation in the repeated i,k in the l.h.s. is over Z). 
The operator Gi contains the creation operators only while G^ contains annihilation 
operators only. Note also that the two operators under the normal ordering commute 
with each other. It is not difficult to find explicitly which rotations are performed by the 
elements Gi,G2,Gz- For Gi^s this is especially simple: 

I n i i'n + Anbi'b , n <0 ^/^n, n <0 

{ ijjn, n>0, { + A^iipi , n > 0. 

For G2 we write G2 = iG^Gg : with G^ = e^-^^-^^ , = e^^s'^^^s. When we move 
through this element, we can ignore either or G2 depending on whether n is negative 
or positive. The rest of the calculation is similar to the one with the normal ordering 
^ (...) ^ given above. It gives: 



G'2^/'„ = (V'n - Anb'^Pb)G2, n > 0. 

It is instructive to write these rotations in the block matrix form: 

i^.G, \^( I A-^\( Gi^. 
) \Q I )\ Gi^+ 

z^_G3 \^( I \ / GsV^- 
^+G3 ) \At I )\ G3^+ 

^_G2 \ ^( I+AZ \( G2^- 

^+G2 j I (/-A+)-l i I G2^H 



(2.36) 
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where the self-explanatory notation is used (we assume that the matrix /— is invert- 





ible). In this notation P+ = 
product of these three: 

/ 



/ 
/ 



/ 



. The full rotation matrix is then obtained as the 



R 



i+A: 







{I-A\ 



I 

At I 



A-_+A-{I-AX)-^At A-{I-AX)-^ 
{I-AX)-^At 



One can check that the second matrix in the last line (which is i? — / = 5) is exactly 
{I—APX)^^A, in agreement with f l2.34p . It remains to calculate the scalar factor in front 
of f l2.33p . Let us take the expectation value of the both sides with respect to the bare 
vacuum. Then we should prove that 

(oo| : exp(^AifcV'r^fc) : |oo) = det(/-P+A). 
Using the decomposition of the operator in the l.h.s., we write: 



E 



{po\ \e 



k>0 



k\ 



Aaibi 



Aa^k^ (00| Vbi • • • ^b.V'l • • • V'ljoo) 



E 

fc>0 



k\ 



E 

a,i,...,a/^>0 



A 


A 


. A 


A 


A 


. A 


A 


A 


. A 



a2ak 



det(/-A+) = det(/-P+A). 



Proposition 12.31 can be easily generalized to the normal ordering *(...)*„. For exam- 
ple, for G given in fl2.27p it can be written as 



G= -exp{Bik^*tP, 



det(/+P>„P) : exp Aik^p*^pk :„ 



(2.37) 



with the matrices A, B determined by the matrix b in (12.270 according to the formulas 

m 

(2.38) 



3 = 6"-!, B-A = AP>nB . 

Here P>„ is the projector on the > n mode space {{P>n)ik 
otherwise). 



6ik for i,k > n and 



2.5.3 Non-invertible (group-like) elements 



The normal ordering allows one to represent in the form (I2.37P not only group elements 
but also some non-invertible elements of the Clifford algebra that satisfy the commutation 
relations 

= Yl ^In^l^ ' i^nG = J2 RnlG^l , (2.39) 
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or 



= J2 R'lnG^i = Yl G (2.40) 



with some (not necessarily invertible) matrices i?, R' (if they are invertible, then both 
pairs of these relations hold with R'R = I but for non- invertible elements only one pair of 
these relations holds). Any normally ordered exponent G = •^e^^-^'^i^k satisfies (12.401) 
and, vice versa, if G is a solution of (12.401) then it can be represented by a normally 
ordered exponent of this form. 

We call elements G of the Clifford algebra such that the commutation relations (12.391) 
or (12.401) hold the group-like elements (in the next section this definition will be further 
extended). If the matrix R fails to be invertible, so does G, as an element of the Clifford 
algebra. In this case it can not be represented in the exponential form (12.271) . Some 
of such elements still admit a representation as a normally ordered exponent of bilinear 
forms in the fermionic operators. However, there are normally ordered exponents such 
that neither R nor R' matrices exist for them. 

Examples. 

a) Let \I', $* be arbitrary linear combinations of the fermion operators ipn) V"^, respec- 
tively, and consider the element 

Q ^ _ , X = 1 + = 1 + a7 - 

where 7:=(oo|\l/<l'*|oo) and a, /3 are related as e^^ = 1 + 7a. For general values 
of a this element is invertible and the two representations are equivalent. However, 
for 77^0ata; = — 1/7 the invertibility breaks down and the element G becomes 

xT/ <T)* 

G 



(oo|^<l>*|oo)' 



which can not be written in the exponential form (I2.27P but can be represented as 
the normally ordered exponent (so that the matrix R exists, but it is not invertible). 
Group-like elements of this type are used in section 13.4.11 below. 

b) For any element G of the form G = I exp (^Aikip^ipk^ ; one can find 

R = I +{I -AP+)-^A and R' = I - {I + AP_)-^A, 

where = I — P+. Then for G = ; exp(^^plipi — ip'^^ip_i^ ; = iplipiip^iip'i^ both 

matrices I—APj^ and I+AP^ are degenerate, so neither R nor R' exist. Nevertheless, 
this element G obeys the BBC (12.411) below. Various projection operators, described 
in section 12. 8[ are also elements of this type. 

2.6 The basic bilinear condition 

It is easy to see that any element satisfying (I2.39P or (12.401) obeys the commutation 
relation 

^kG ®rkG = J2 ^"^k ® Grk (2.41) 
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which we call the basic bilinear condition (BBC). It means that G ® G commutes with 
Tlk'^k ® "ipk- terms of matrix elements it reads 

J2 {U\ ^kG \V) {U'\ r^G \V') = J2 {U\ G^, \V) {U'\ Gr^ \V') (2.42) 

for any states \V) , \ V'), {U\ , {U'\ from the space T-Lp and its dual. Indeed, substituting 
f l2.39p or f l2.40p instead of i/jkG and Gi/j^ or i/j^G and Gipk in the left and right hand sides 
of fl2.4ip respectively, we get the identity. 



Proposition 2.4 All normally ordered exponents of bilinear forms solve the BBC. 

It is enough to prove the statement for the normal ordering with respect to one particular 
vacuum, say |0). Let us consider the decomposition fl2.35p of general normal ordered 
exponent. Then, Gi and G3 are group elements, so they satisfy BBC and all we have to 
prove is that G2 also does. This is obvious from fl2.36p and corresponding rotation for 
rk's: 

( i)nG2 = G2{iJn + ^nSV^fe), ^ < ( Ga^^ = {i^* + A^iijl)G2, n<0 

[ GaVn = ii^n - A^f,iJb)G2, n > 0, \ ^^Ga = G2(C - A^bi^l), n>0. 

(2.43) 

Indeed, one can easily check that ( I2.4ip is valid separately for the sums with A; < and 
k>0. 

It turns out that in addition to the group elements and normally ordered exponents 
there are other solutions to the BBC. For example, it is easy to check that G = as 
well as G = ^Z'* solve f l2.4ip . More generally, so does any linear combination of iJj^s (as 
well as '?/'*'s). Indeed, set v = J^k'^k'ipk, then 

^ vipk '^vipi = - ^ i)kv ® {-i^lv + Vk) = ^ i'kv '^i'lv Vki^kV ® 1. 

k k k k 

The last term vanishes because Y2k "^ki^kV = f ^ = 0. At the same time, the element 
G = ipn (for example) does not generate any linear transformation in the linear space 
spanned by ipl, neither of the form fl2.39p nor fl2.40p . 

It appears that the only important property of an element G for what follows is not 
the induced rotation like fl2.39p or f l2.40p but the BBC fl2.4ip . Therefore, from now on 
we extend the definition of the group-like elements including in this class all solutions to 
the BBC. 

We have the following general properties of the group-like elements. 



Proposition 2.5 The elements that satisfy the BBC ^2.4-1^ form a semigroup: if G and 
G' satisfy it then so does GG' . 

This is obvious from fl2.4ip . 
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Proposition 2.6 All solutions of the BBC ^2.41^ have definite charge, i.e., [Q,G] = qC 
with some integer q. 



The proof is given in Appendix B. 



2.7 The generalized Wick theorem 



Let Vi = '^jVijipj be linear combinations of ip/s only and w* = '^jW*j'ip* be linear 
combinations of ip^s only, as before. 



Theorem 2.7 Let G,G' be any two group-like elements with zero charge. Then for any 
Vj, w* and any n such that {n\ G'G \n) ^ it holds 



{n\ 


G'vi ... 




. . . wlG 


\n) 


{n\ 


\G'G\ 


\n) 





G'vjW*G 


\n) 




\G'G\ 


\n) 



det ^'^ (2.44) 



This is the generalized Wick theorem. A similar statement holds when w*^s stand to 
the left of v/s. Writing {n — m\ = {n\ ipn-i ■ ■ ■ i^n-m or {n + m\ = {n\ ip^ ■ ■ ■ "^n+m-i with 
m > 0, we get from (I2.44p : 

Corollary 2.8 Let G be any group-like element. Then 

..det W;~MG\n) ^ ^2.45) 



{n — m\ 


w 


,* 

m ■ 


..wlG\ 


\n) 


{n\ 


\G\ 


\n) 



{n + m 


V 


'm • 


..viG\ 


\n) 


{n\ 


\G\ 


\n) 



(n 




ViG 


\n) 


{n\ 


\G\ 


\n) 



det ' ' ^."7" , . ' ' . (2.46) 



The theorem can be proved by induction using the BBC in the form f l2.42p . Suppose 
f l2.44p is valid for some m > 1 (it is trivially valid at m = 1). Set 

{U\ = {n\G'wl, {U' \ = {n\G' viV2. . .Vm+iw*^^]^uj*^. . .wl, \V) = \V') = \n) . 

Plugging this in the BBC, we see that its r.h.s. vanishes because either ipk \n) = or 
= while the l.h.s. gives 

^ {n\ G'wlipkG \n) {n\ G'vi . . . Vm+iw*^^^ . . . wli)lG \n) = 0. 

k 

Substituting wl^pk = — i'k'wl in the first factor and moving ipk the second factor 
through the chain of the w*^s, we get in the l.h.s. 

{n\ G'G \n) {n\ G'vi . . . Vm+iw^_^^i . . . wlG \n) 
- i-ir Yl (^1 G'^kwlG \n) {n\ G'v, . . . Vm+irkW*m+i ■ ■ ■ \n) . 
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We proceed by moving ijjl to the left through the chain of the Wj's using Vjipl = Vjk — '4^k^j 
at each step. The result is 

{n\ G'G \n) {n\ G'vi . . . 'Um+i'W^m+i • • • ""^i^ l'^) 

m+l 

'vjwlG \n) {n\ G'vi ... . . . fm+iu^m+i • • • 1^) 

i=i 

+ ^ {n\ G'tpkwlG \n) {n\ G'iplvi . . . Vm+iw^_^-^ . . . w^G \n). 

k 

In the last line we can again use the BBC to write it as 

(^1 i^kG'wlG \n) {n\ 'iplG'vi . . . Vm+iwl^^^ . . . w^G \n) 

k 

which is because either {n\ ^/'fc = or {n\ tpl = 0. Therefore, we conclude that 

{n\ G'G \n) {n\ G'vi . . . Vm+iw^_^^i . . . wlG \n) 

m+l 

+ ^^(—1)-' {n\ G'vjwlG \n) {n\ G'vi ... ffj ■ ■ ■ fm+iWm+i . . . W2G \n) = 

i=i 

or 



{n\ 


G'vi . . . v„,+iw*^^^ . . . wlG 


\n) 


{n\ 


\G'G\ 


\n) 



{n\ 


G'vjWlG 


\n) {n\ 


G'v^ . . . ^, . . , 




<+i • • • 


\n) 


{n\ 


\G'G\ 


\n) 






\G'G\ 


\n) 



{n\ 


G'vi . . . • 


VmG"w*^ 


. . . wlG 


\n) 


{n\ 


G'G"G 


\n) 



m+l 

By the assumption, the second ratio in the r.h.s. is the m by m determinant. We see 
that the r.h.s. is the expansion of the m + 1 by m + 1 determinant in the first column, 
so the theorem is proved. 

In a similar way one can prove that for arbitrary group-like elements G, G', G" it holds 

{n\G'v^G"w*G\n) 

Moreover, since for the derivation only the BBC (12.411) is used, this version of the Wick 
theorem can be immediately extended to solutions with non-zero charge. Namely, for 
any three solutions of (12.411) Gq, G'^, and G'^„ with charges q, q' and q" respectively, one 

ln\ G' .lu . . . i).^.G"..uf . . . ui*G. I«\ 

det 

i,jr'=l,...,rri 

where n = n — q — q' — q" . For a general version of the Wick theorem see also [H]. 

There exists an alternative determinant representation of the expectation value in the 
l.h.s. of (I2.45p . which is another form of the generalized Wick theorem. 





G',,v, . . , 


■ nmG'^„w;^,. 


. . wlGq 


\n) 






\n) 



{n\ 


G'giVjG'^„W*Gq 


\n) 


{n\ 


f^i (-ill 
^q'^q"^q 


\n) 
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in, — 777,1 in* in^Cr \n) 




(77 G 77) 




(77 + 777 f m • • • viG \n) 




(n G 77) 


= 1 


(77 Gvi . . .Vm\n — m) 




(77 G 77) 


= I 


{n \ Gwl ■ ■ ■ \^ + 




(77 G 77) 





Corollary 2.9 Let G &e any group-like element. Then 

det + (2.48) 

:=l,...,m (77-J + ll G I77-J + I) ^ ^ 

det (n+j\-^G\n+j-l) ^2.49) 
=i,...,r?i (n+j — 1| G \n+j — l) 

=l,...,m (77— J + l| G |77— J + 1) 

det (-+J-^\Gn^nn+j) 

:=l,...,m (77+J-ll G |77+j-l) ^ ^ 

Let us outline the proof of fl2.48p . The idea of the proof is as follow^]. The first 
columns of the matrices in the r.h.s. of fl2.45p and fl2.48p coincide and one can show that 
the j-th column of the matrix in fl2.45p is equal to the j-th column of the matrix in fl2.48p 
plus a linear combination of the first j — 1 columns. To show this, we use the BBC. Let 
us take 

{U\ = {n + l\iPiw*, {U'\ = {n\, \V) = \n), \V') = \n + 1) (2.52) 
in fl2.42p . It is easy to see that only one term in the r.h.s. remains: 

^ (77 + 1| i/jiw-ipkG \n) {n\ iplG [77 + 1) = (77 + 1| ipiw*G \n + 1) (77] G I77) . 

k 

In the l.h.s. we plug w* = w*jip* and move ipk in the first factor to the left position. 
Using the anti-commutaton relation for the fermion operators, and transforming the sum 
J2k '^iki'k back to w* in the second factor, we arrive at the 3-term identity 

(77I G |n) (77 + 1| -ipiw^G |77 + 1) = (77 + 1| G |n + 1) (n| ■ipiw*G \n) 

+ (77 + 1| ipiG \n) {n\ w*G |77 + 1) . 

Divide it by (77I G [77) (77 + 1| G |77 + 1), take a sum over the values of 77 equal to 77 — j, 77 — 
j + ..,77 — 1, and then put I = n — j. We obtain the relation 

{n\'4)n-j'w*G\n) _ {n-j\w*G\n-j + l) {-X {n-k +l\ilJn-jG\n-k) {n-k\w*G\n-k+l) 
{n\G\n) ~ (n-j +1| G|n-j +1) ^ {n-k\G\n-k) {n-k+l\G\n-k+l) 

which shows that the columns of the matrix in f l2.45p are the linear combinations of 
the columns of the matrix in fl2.48p and their determinants are thus equal. The proof of 
f l2.49p is similar. 

Again, formulas fl2.48p - fl2.5ip can be easily generalized to the case of group-like ele- 
ments with non-zero charge. 



^This proof was suggested by Z.Tsuboi. 
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2.8 Projection operators 

In some applications (in particular, to models of random matrices), an important role is 
played by the non-invertible group-like elements 

p+ = :exp ( ^v^.v^n : = Hii-r.i^,) = HV'^^r^ 

\i<0 / i<0 i<0 



P' = : exp (-J2^^i^n : = 11(1 - ^,r^) = H v^*^- 

V i>0 / i>0 i>0 



(2.53) 



In a sense, these operators are projectors to the spaces of states with positive (i.e., 
> 0) and negative (i.e., < 0) charge respectively. Their properties (extensively used in 
what follows) can be easily seen from the definition. Both operators obey the projector 
property: (P^)^ = P^. The operator P+ kills negative creation modes standing to the 
right and negative annihilation modes standing to the left and commutes with all positive 
modes: 

P+i;l=i^kP^ = 0, k<0, 

(2.54) 

[P+,^*] = [P+,^,]=0, k>0. 

The operator P~ kills positive creation modes standing to the right and positive annihi- 
lation modes standing to the left and commutes with all negative modes: 

p-^k=^lP~ = 0, k>0, 

(2.55) 

[p-,rk\ = [p-,^k] = o, k<o. 

From this it is obvious that P^\n) = at n < and P'''|n) = \n) at ri > 0. Similarly, 
P-|n) = at n > and P~\n) = \n) at n < 0, so that P+|0) = P-|0) = |0). Somewhat 
less obvious properties, also used in what follows, are 

p+e-^-W |n) = ^ (-l)''WsA(t)|A,n), n > 0, 

e{X)<n 

(2.56) 

(n|e-^+(*)p+ = J] (-l)^WsA(t) (A,n|, n > 0, 

l{X)<n 

(see section [2l9] for definitions) in particular, P+c^" |0) = |0), (0| c^+P^ = (0|- 

Extending definition f l2.53p to other vacuum states, we can also introduce the projec- 
tors 

i<n i>n 

with similar properties. Moreover, similar operators can be introduced for any basis 
state: 

oo 

k=l k&In,\ 

(2.58) 

oo 

Pn,\ = X\'^*n-\'^+k-l'^n-\'^+k-l= H ^^^^ 
k=l HI„_x 
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(in particular, P+ = P^0). They also obey the projector property (P^^)^ ~ 
commute with each other, [P^A'^nA] ~ 0. Properties fl2.54p and fl2.55p are naturally 
generalized for the corresponding creation and annihilation operators. Their product 
^nx^nx projector to the state \n, A): 



(2.59) 



Thus for the basis states we can write |A, n) (A, n\ = P^x^n x ^^^^ particular, |?7,) (?t,| = 
P+P- and |0) (0| = P+P-. Note that 

oo 

Y^\x,n){x,n\ = Y,p-^,pi,= n {i^krk + rM = l, 

n,X n,X k=—oo 

as it should be. 

From fl2.16p and (12.170 it follows that the operator \n, A) {m,fi\ is represented by an 
infinite product 



|n. A) (m,/i| = Yl Rk 



(n,A;m,^) 



(2.60) 



fc=— oo 



where R 



(n,A; m,/i) 



is given by 



(n,A; m,/^) 



' iJki^i, k e in,x, k G /, 

i^llpk, k^In.X-, k^ Im,,, 
Ipk, k G /n,A) k ^ Im,fj, 
^ ^fc' ^ ^ -^n,A; k G Im,fi ■ 

Any operator \n, A) (m, /i| is group-like. 



(2.61) 



2.9 Currents J±(t) 

Among bilinear combinations of fermions of particular importance are the operators 

Jk = Yl I'P^'pHk: = res,(:^(z)^^-v*W:) (2.62) 

which are Fourier modes of the "current operator" J{z) = ',ip{z)il!*{z)l. Note that the 
normal ordering in (I2.62p is essential at = only, and in this case Jq coincides with 
the charge operator Q (12.241) . At A; 7^ the normal ordering can be ignored: 

These operators are of the form (12. 3p with the matrix Aij = Sij-k- This matrix is not 
of the class discussed in section 2.1 because it has an infinite non-zero diagonal, and 
thus one should take some care when working with formal expressions containing infinite 
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sums (see the example below in this section). It is easy to see that [Jk,'ipm] = 4'm-k-i 

Let tfe, /c e Z, be arbitrary parameters (called times). It is convenient to denote 
the collection of times with positive (negative) indices by t+ = {^i,^2,---} and t_ = 
^_2, . . .} respectively. Set 

J+ = J+(t+) = ^4Jfc, J_ = J_(t_) = J]i_feJ_fe (2.64) 

k>l k>\ 

and introduce the generating function 

k>l 

Here and below we write simply t for any half- infinite set of times (either t+ or t_) when 
it does not lead to a misunderstanding. 

It is easy to check that the series ip{z), ip*{z) transform diagonally under the adjoint 
action of the elements e'^+, e^^: 



e'^+(*V(^)e"^+^*^ = e«(*'^V('2), 
gj+(t)^*(^^)g-j+(t) ^ e-^(*'^V*('2) 



(2.65) 



and 

e"^-W^(^)e--^-W = e^(*'^"')V('2), 



(2.66) 



gj_(t)^*(^)g-j_(t) ^ e-^(*'^"')V*('2)- 
In terms of the polynomials hk{t) defined by 

e«(*'^) =^/ife(t)^'= (2.67) 

fe>0 

the adjoint action of e''* to V'n, -0* is given by the formulas 



fe>0 fe>0 



(2.68) 



fc>0 fc>0 

Their easy consequence is 

(l|e-^+W^,|0) = /i,(t), j>0, (2.69) 

and, more generally, 

(0|V.V+(%-|0) = V^(t), hj>0- (2-70) 
The commutator [J^, J;] can be calculated as 

j 3 
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Using the formulas obtained above, we get 

[Jk. Ji] = i^i-i^^j+i - MUi+k) ■ (2.71) 
j 

Naively, one can shift the summation index j — )■ j + fc in the first sum to get 0. In fact this 
can be done only when + / 7^ 0. At k = —I the r.h.s. appears to be ill-defined because 
of the infinite summation. Let us fix the r.h.s. of equation (12.711) . When + / 7^ 0, the 
r.h.s. is well-defined (and equal to zero) as it stands because the result of its action to 
any basis state contains only a finite number of terms. When k + I = 0, one should first 
rewrite the r.h.s. in terms of the normally ordered expressions: 

where 6{j < k) = 1 for j < k and otherwise. The normally ordered expressions are 
well-defined and the summation index can be shifted thus yielding 0. The rest gives the 
commutation law 

[Jfc, Ji] = k5k+i,o ■ (2.72) 

We see that the Fourier modes of the fermionic currents commute as bosonic operators. 
We thus have 

[j+(t+), j_(t_)] = j2ktkt-k 

k>l 

and 

gj+(t+)gj_(t_) ^ exp(j2ktkt-k)e'-^'-^e'+^'+\ (2.73) 

k>l 

2.10 Expansion of the states e^*^*^ |A, n) and Schur functions 

The definition of the vacua implies that {n\ J_ = 7+ \n) = and thus {n\e'^~ = {n\, 
6"^+ \n) = \n). The coherent states e'^~ \n) and {n\e'^+ (and, more generally, the states 
\X,n)) can be expanded [2l |12] as linear combinations of the basis states. The coeffi- 
cients are the famous Schur polynomials. This expansion is important since it provides a 
link between hierarchies of integrable equations and the theory of symmetric functions. 

2.10.1 Expansion of the states e"^**^*^ \n) 

Given a Young diagram A = (a|/3), one can introduce the Schur polynomials (or Schur 
functions) via the Jacobi-Trudi formulas [9]: 

s,(t)= det /iA.-.+,(t), (2.74) 

«J = lv,^(A) 

where hk{t) are defined in (12.671) (they are Schur polynomials for one- row diagrams). For 
details see [9] and Appendix A. We recall that 6(A) = + 1) (see (12151) 1 For the 

empty diagram S0(t) = 1, 6(0) = 0. 
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Proposition 2.10 It holds 

e^-(*)|n)=5^(-lp)s,(t)|A,n), (2.75) 

A 

(n|e^+(*) = J](-l)^Ws,(t)(A,n|, (2.76) 

A 

where the sums run over all Young diagrams A including the empty one. 

From orthogonality of the basis states it follows that fl2.75p and f l2.76p are equivalent to 

(A,n|e^-(*)|n) = (-l)^WsA(t), (2.77) 

(n| e-^+(*) |A, n) = {-if^^^ s x{t) . (2.78) 

Let us sketch a proof of the first formula of the proposition. First of all we recall that 
the Schur polynomials for an arbitrary Young diagram A = (a|/9) are expressed through 
the Schur polynomials S(o|^)(t) corresponding to the hook diagrams A = (a + 1,1'^) with 
the help of the Giambelli formula 

SA(t) = .. det s(«.,/3,)(t). (2.79) 

i,j=l,...,d(\) 

Therefore, the r.h.s. of fl2.75p can be transformed as follows: 
5^(-l)''W..(t)|A,n) 

A 

= Yl Yl (-1)''^^^ l<1^-</("^l'^^)^*^ V^n-/3i-l • • • ^n-/3,-l^n+a, • • • ^n+aj^) 

d>0 ai>ti2> - >"d>0 
/3i>/32>...>/3jj>0 

/3l,/32,...,/3d>0 



exp(^ ^ (-l)^^^S(a|/3)(tX_^_lV'n+a) 1^) ■ 



a,l3>0 

Next we use the relation 

S(a|/3)(t) = (-1)^ Y ^/3-m(-t)/ia+m+l(t) 
m=0 

which is a consequence of the Jacobi-Trudi formulas (see Appendix A). With the help of 
equations (12.680 we can write, continuing the chain of equalities above: 



5^(-l)^W.,(t) |A, n) = exp(- Y rM-it)) {^rn{J-m>„) \n) 
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where 

fc>0 

(2.80) 

fc>0 

and (V^m(</-(t)))>„ means that in the series ipm{J~i't)) only T/^-operators ipk with k > n 
are kept. The next step is to notice that the operator in the r.h.s. has the same effect 
when acting to the right vacuum as the normally ordered operator 



(which was considered in so 

5^(-l)^W..(t) |A,n) = ^ exp(- ^C(^-(t))^™(J-(t))) ^ In) 



m<n 



e^-(t) X exp(- Yl "P^mi^rn) e'^'^*) \n). 



(A comment on the last step is in order. The substitution -ipmiJ-it)) = e'^"^*Vme~"^~^*\ 
ipl^{J-{t)) = e'^-'^^Vm^"'^"*'*'' be done under the sign of normal ordering. It can be 

done only after expanding the normally ordered exponent in a series of the form f l2.25p .) 
Finally, it is easy to see that the operator 



exp(-^ ^ = ^ n X = ^ n (1-^-^' 

m<n m<n m<n 

(2.81) 

n(i-c^rn)= n^™^™ = 



(see fl2.57p ) has the property P^e \n) = \n), which proves the desired formula. 
Using the Cauchy-Littlewood identity ( \A7\i . we can write 



e-'-^*) |n) = exp J^ktk) \n) = ^ sx{t)sx{3-) \n) , 

k>l A 

hence 

\X,n) = i-lf'^sxi3-)\n) (2.82) 

and, similarly, 

(A,n| = (-l)''W(n|sA(J+), (2.83) 

where J± = (J±i, ^±2/2, ^±3/3, . . .). These formulas for the basis vectors were obtained 
in [13]. 
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Because of importance of Proposition 12. 101 we find it instructive to give another proof. 
Instead of proving fl2.75p we will prove the equivalent formula f l2.77p : 

(A, n| e-^-^*) \n) = {n\ . . . V^n-/3,(A)-i • • • i^n-h-i^^'^'^^^ I"-) 

= det {n\ rn+a.i^n-p,-ie'-^'^ \n) 

i,j = l,...,d(\) 



det (n|C+a,(^-(-t))V^n-/3,-i(J-(-t))|n) 

i,j = l,...,d(\) 



(2.84) 



det V/ic,,-fe(t)/i/3^+i+fc(-t) 

ij=l,...,d(A) 
= (-l)^Ws,(t), 

where we have used the Wick theorem f l2.44p . the expression for the Schur polynomial 
corresponding to a hook diagram (lASP and the Giambelli formula fl2.79p . 

2.10.2 Expansion of the states e"^**^*) in the skew Schur functions 

The previous proposition can be generalized as follows: 

Proposition 2.11 It holds 



e-^+(*)|A,n) = 5^(-l)^(^\'^)s,v.(t)|^,n), 
(A,n|e-^+W = 5^(-l)''('^\^).,\,(t) (^,n|, 
(A,n|e''-(*) = 5^(-l)''(^\^).Av(t) (/x,n|, 



(2.85) 



where 6(A \ /i) = 6(A) — 6(/i) and 

SAv(t) = det /iA,-M,-.+,(t) (2.86) 

are skew Schur functions (see and Appendix A). 

Note that the sums in the second and the fourth lines are finite because SA\^(t) is non-zero 
only if C A. These formulas imply 

(/x,n|e^+(*) \X,n) = (A,n| e^-^ \fi,n) = {-lf^'\^h^\,{t). (2.87) 
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Let us sketch the proof of (12.851) . The first formula is proved by the following chain of 
equalities: 

e'^-(t)|A,n) = (-l)''We^fc>i*'=^-*SA(J-) In) 

= (-l)^W5^s.(t)s,(J_)sA(J_)|n) 

V 

where in the second line the Cauchy-Littlewood identity is used. In the third line we use 
the definition 

V 

where the Littlewood-Richardson coefficients c^^ are determined by 

S^(t)s.(t) = 5^c;,SA(t). 

A 

The other formulas are proved in a similar way. 

The determinant formula (12.861) follows from the Wick theorem. Put I = i{X) and 
add £(A) — £(/i) zero lines to the bottom of the diagram if i{fi) < i{X). Using (I2.13p . 
we write: 

SAv(t) = r,,-i . . . r,,-ie'^^'^^x,-i . . . ^x,-i 

= det {0\re+,,-,e'■^^'^^Pe+x.-^\0) 

= det /iA,-^,-i+j(t). 

l<i,j<l ^ ^jy ^ 

The last equality follows from (I2.70p . 

2.11 The boson- fermion correspondence 
2.11.1 Bosonization rules 

Consider the current operator 

As we have seen, its Fourier modes Jk have the same commutation relations as bosonic 
operator modes: [Jk, Ji] = k6k+ifi, so they generate the Weyl algebra and Jk and J-k/k 
are canonically conjugate. 
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The operator Jo = Q is special. Its canonically conjugate partner is an operator P 
such that is the shift operator acting as 

p I -p I p I * -p I * 

on the fermionic operators or e"^^ \n) = \n ± 1), {n\ e^^ = (n =p 1| on the vacuum states 
(this definition imphes the commutation relation [Q, P] = 1). In accordance with Propo- 
sition \2.2\ 

n,X 

where the states in r.h.s. are given by (ICTll . Note that e^^ satisfy the BBC. 
All this prompts to introduce the chiral bosonic field 

k>o k>o (2.88) 

= J4[z]) + P + Jolog z-J+{[z~']) 
such that zdz(f){z) = J{z) or 

/■""^ (]z 
<f>{z,) - <f>{z,) = / J{z)-. 

In the second line of fl2.88p we use the notation ^±([2:]) = J±iz + ^J±2Z^ + | J±3-2^ + . . .. 

The operators J+fc with A; > kill the right vacuum, so they are bosonic annihilation 
operators while J-k with A; > kill the left vacuum, so they are creation operators. 
Let us introduce the bosonic normal ordering *(...)* defined by the usual rule that all 
annihilation operators are moved to the right and all creation ones to the left and consider 
the normally ordered exponents of the free bosonic field: 

(2.89) 

(note that the normal ordering acts only to the current modes Jk with nonzero k). From 
f l2T3|) it follows that 

gj±{t) ,^4>iz)* g~J±(t) = gC{t,^±i).g0(^). 

gj±{t) ,^-4>iz)* g-J±{t) ^ g-5(t,2±l).g-</.{2). _ 

These formulas tell us that •e^''^^^-* * behave like the fermionic fields 'ip(z),^*{z). Moreover, 
one can show that all matrix elements of *e^'^*^^) * with respect to the basis states coincide 
with those of iIj{z),iIj*{z). Therefore, we can identify 

^(z) = :e^(^):, ^*(^) = :e-^w:. (2.90) 

Let us outline the proof of the first formula (the second one is proved in a similar way) . 
We need to compare matrix elements of the operators in both sides. It is clear that the 
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matrix elements of ip{z) vanish unless the charge difference between the left and right 
states is 1. Consider the generating function of non-zero matrix elements, which is the 
matrix element between the coherent states {n\e'^+^^+^ and e*'"^*-^ \n), then we should 
prove that 

^^|gj+(t+).g<^(.),gj_(t_) \n-l) = (n|e*^+(*+V(^)e"^-^*-^ \n - 1) 

for all t-i- and n. It is a matter of straightforward calculation, which uses the explicit 
form of le'^^^'^l (12.891) . the commutation relation (12.731) and relations (I2.65p . (I2.66p . to 
see that the both sides are equal to 

z^-' exp (e(t+, z) - e(t_, z-') + J2 ktkt-k) ■ 

k>l 



Let us check that the bosonization formulas (I2.90p imply the relation 

'Mz)riz): = zd,(i)iz) (2.91) 

and in this sense are consistent with (12.881) . We start with ip{z2)ip*{zi) = ♦e'^^^^^ * *e~^^'^'^^ * 
and rewrite the both sides through normally ordered expressions. An easy calculation 
shows that 

^(Z2)^*(^i) = 'Mz2)rizi): + ^^, 

Z2 - Zi 

.g</'{22)*,g-</'(2l), ^ ,g0(z2)-</'(2:i),_ 

Z2- Zi ' 

Setting Zi = z, Z2 = Zi + e, we thus get: 

-l) = i^fj{z + e)^lj*{z): 
which coincides with (I2.9ip in the limit e — )■ 0. Next terms of expansion are given by: 

'Mz + e)riz): = z(^d<p + ^{: (s^)' : + sv) 

(2.92) 

+ |y (: (50)' : + 3:909^0: + 9=^0) + o(e=^) j . 

where d"^(t> := d'^chiz). For more details see [1 



Being applied to vacuum states, the operators (I2.89P simplify because either e^"^"^'^'^ 
or e^'^+*^[^ l-* disappears. Acting by both sides of (I2.90p to the left vacuum, we get the 
bosonization rules 

(2.93) 

(n| ^*{z) = 2-" + e^+([^"'l). 

Their simple consequence is 

{n\ ^*(C)^(2) = — ^ {n\ e-^+([^"'l-[^"'l). (2.94) 
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The repeated use of equations (12.931) gives the following formulas for arbitrary number 
of fermionic operators: 



(n| . . . = (zi . . . z,r.)-"^^^{z.--Zj) (n+m| e^+(I-r'l)+-+^+(fc'l). 



(2.95) 



(2.96) 



Merging the points in these formulas, one gets 

(n| a™~V(^) • • • d'^{z)i){z) = a^z™("-") {n - m\ e~"-^+([^"']), 

{n\ d"'-^ij*{z) . . . dij*{z)ij*{z) = arn^"""^""^™"^) {n + m\ e'^-^+d^"'])^ 
where = 1! 2! . . . (m — 1)! 

Similarly, the bosonization formulas for the right vacuum are: 

ipiz) \n) = ;z"e-^-(W) \n + l) , 

ij*{z) \n) = ;2-"+ie--^-(W) 
from which it follows that 

^*{C)tp{z) \n) = e^-(W"[^]) \n) (2.98) 

and 

tP{z^) . ..ij{zj \n) = {z, . ..z^Y \{{z, - z,) e^-([-il)+-+^-(M) |n+m). 



(2.97) 



r{z,) . . . r{zm) \n) = {z,... ;^„)-"+i l[{z, - z,) e-^-([-i])-...-J-([^™]) |^_^)^ 



(2.99) 



i<j 



ip{z)dip{z) . . . (9'"-V(^) 1^) = a^2'"("+™-^)e"^'^-(W) |n + m) , 

(2.100) 

^p*{z)d'ip*{z) . . . 9™-V*(-2) I"-) = a^z'"(™-")e-'"-^-(W) \n-m) . 
2.11.2 Vertex operators 

Let us present the bosonization formulas in a different but equivalent form which is based 
on an explicit realization of the bosonic Fock space "Hb by polynomials in infinite number 
of variables ti, ^2, ^3, • • •• More precisely, following [3], consider the space 

Kb = C[w, w-\ h, t2, t3, . . .] = C[ti, t2, t3, . . .], 

where an extra variable w is added to take into account the fermionic states with different 
charges, and the map ^ : T-Lf ^ 'Hb defined for an arbitrary state |f/) G "Hf as follows: 

$(|f/)) = ^y;'(/|e-'+(*) (2.101) 
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If the state \U) has a definite charge g, then the sum in the r.h.s. contains the only non- 
zero term with I = q. It can be shown that this correspondence means that the fermionic 
Fock space Hp is isomorphic to the bosonic Fock space T-Lb and the Fourier components 
of the current as well as fermionic creation and annihilation operators become some 
operators acting in the space of functions of w and tj. 

Proposition 2.12 For the components of the current we have: 



wd^<i>{\U)) , k = 0, (2.102) 
-H.^mU)), k<0. 



The first two formulas are obvious from the definition f l2.10ip . The third one follows from 
the commutation relation 6'^+*^*^ J_fc = (J_fc + ktk)e'^+^^^ (where k > 0) obtained with the 
help of formula (12. 4p . 

Let us introduce vertex operators by the formulas 



X(z) =exp [J2t,z^j exp (^-J^J^^^/j ^''''^ 



X*{z) = exp I - 5^ j exp ^*.) ''''' 



(2.103) 



where the operators P, Q are defined by their action to the functions of w: 

e^fiw) = wfiw) , z^fiw) = f{zw) 

(it is easy to check that [Q, -P] = 1). Using the previously introduced short-hand notation, 
we can write: 

X{z) = e«(*'^)e-«(^'^")e^z«, 

(2.104) 

where d = {dt,, ^dt^, |(9t3, . . .}. 

Proposition 2.13 Under the boson-fermion correspondence, the fermionic operators 
ip{z), 4'*{z) are represented by the vertex operators X{z),X*{z): 

^^Piz)\U))=Xiz)^\U)), 

(2.105) 

^r{z)\u)) = x*{zmu)). 

In order to prove the first equality, we write down the both sides separately, 

$(^(^) \U)) = J^m;" (n| e^+(*V(^) \U) = e«(*'^) ^m;" (n| ^{z)e^+^''> \U), 

n n 

X(z)<l>(|f/)) =e«(*'")^;z"-iw7"(n-l|e"-^+([""'])e^+(*) |f/), 

n 

and apply the bosonization rule (I2.93p . The second equality is proven in a similar way. 



32 



3 Tau- functions 



As it has been established in the works of the Kyoto school, the expectation values 
of group-like element^ are r-functions of integrable hierarchies of nonlinear differential 
equations. This means that they obey an infinite set of Hirota bilinear equations [151 US] • 

3.1 The main types of r-functions 

There are three main types of the r-functions (which generalize each other): 

• The r-function of the KP hierarchy depending on the times t = {ti,t2, . . .}: 

r(t) = (Ole-'+^^^GlO). (3.1) 

• The r-function of the modified KP (MKP) hierarchy ^ [TH [13 [2ojl. 

r„,(t) = {n\ e-^+(*)G \n) . (3.2) 

Equations of the MKP hierarchy are differential-difference equations which include 
shifts of the variable n = to (the "zero time"). At each fixed n, the r-function (13. 2p 
is a r-function of the KP hierarchy, so the ra-evolution represents an infinite chain 
of Backlund transformations. 

• The r-function of the 2D Toda lattice (2DTL) hierarchy [2T| [22]: 

r„(t+, t_) = {n\ e'^+(t+)Ge--^-(*-) \n) . (3.3) 

It is the most general r-function associated with the one-component fermions. At 
fixed n, t_ it is a r-function of the KP hierarchy (as a function of t_|_) and at fixed 
n, t+ it is a r-function of the KP hierarchy (as a function of t_). 

The meaning of the variable n depends on the class of solutions. For some classes 
of solutions n can be regarded as a continuous, or even complex, variable. Using the 
operator P introduced in section I2.11[ we can equivalently represent the n-dynamics as 

r„(t+, t_) = (0| e-^+^^+^e-^^Ge^^e--^-^*-) |0) . (3.4) 

Remark. More general r-functions can be obtained as matrix elements of group-like 
elements with non-zero charge between different Dirac vacua. Any group-like element Gq 
with charge q gives rise to the r-function 

r„(t+, t_) = (n| e-^+^^+^G.e-''-^*-) \n - q) (3.5) 

which, along with (13. 3p . is a solution to the bihnear identity, see (I3.17P below. In section 
13.4.11 we show that a family of (quasi)polynomial and soliton KP and MKP r-functions 
corresponding to group-like elements with non-zero charge can be understood in terms 
of singular limits of group-like elements with zero charge. 

^In this section we assume that the group-hke elements G have zero charge, if the contrary is not 
stated exphcitly. 

■^In a shghtly more general version called there a coupled modified KP (cMKP) hierarchy, was 
considered. 
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3.2 Bilinear equations for the r-function 

3.2.1 The bilinear identity for the MKP and KP hierarchies 

Consider first the MKP r-function. It obeys a bilinear identity which is a direct conse- 
quence of the BBC in the form f l2.42p . Setting \ V) = \n), \V') = \n') with n > n', where 
\n) and \n') are two shifted Dirac vacua, we have 

J^iUlipkGln) {U'liplGln) = (3.6) 

fcez 

because either tp^^ or tpl kills the state in each term in the r.h.s. of f l2.42p . Now, setting 
(f/| = {n + l\ e-^+(*), {U'\ = {n' - 1| e-^+(*'), we can write 

= ^(r2 + l|e'^+(*VfcG'|n)(n'-l|e'^+(*'VfcG'|n') 

k 

= les, \z-^ {n + 1| e-^+(*V(^)G' \n) {n' - 1| e-^+(*'V*(-2)G' 1^') 



res. 



where we have used the commutation relations of the fermion operators with 6"^+*^*^ and 
the bosonization rules (12.931) . Here lesz means picking the coefficient in front of of 
the Laurent series. In this way we arrive at the bilinear identity 



e 



C(t-t',.)^n-„'^^(^ - [z^'jWif + [z~'])dz = (3.7) 



Co 



for the r-function valid for all t,t' and n > n'. Here we use the standard short-hand 
notation 

t±[^] = {ti±z,h±-z'',ts±^z^...}. 
At n = n' (13. 7p becomes the bilinear identity for the r-function of the KP hierarchy. 



e 



r(t - [z'^])T{t' + [z~^])dz = 0. (3.8) 



Co. 



Note that it implies that if r(t) is a KP r-function, then so is r(— t). At n = n' + 1 (13.71) 
gives the bilinear identity for the MKP hierarchy: 

ze«(*-*' ^)r„+i(t - ri])r„(t' + [z-'])dz = 0. (3.9) 

Coo 

The same bilinear identity can be obtained for the r-functions constracted as matrix 
elements of group-like elements with non-zero charge. 

The choice of the integration contour Coo depends on the type of the time evolution. 
Formally, the evolution factor e^''*^* ,z)^n~n -^^^^ essential singularity at oo, and the 
contour should then be a small circle around oo (a big circle in the complex plane). 
This standard prescription does work if only a finite number of times are nonzero. 
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However, in general, when the values of the times are such that the factor e^'^*"*''^)^;"'""'' 
has singularities at finite points of the complex plane, then the prescription is as follows: 
the contour Coo must encircle all singularities of the function r„(t — [2;^^])x„/(t' + [-2"^]), 
leaving all the singularities of e^'-*"* ,z)^n-n outside the contour. 

Remark. There is a freedom to multiply the r-function of the MKP hierarchy by 
an exponent of any linear form of times with constant coefficients and by an arbitrary 
function of n: 

Tn{t) ^ C(n)exp(^Cfetfc)r„(t). (3.10) 



fc>i 



For the KP hierarchy 



r(t)^Cexp(5^Cfet,)r(t). 



(3.11) 



k>l 



Clearly, these transformations preserve the form of the bilinear identities. By noting that 



C7(0):exp (j^rf,^;^, ) : l^) = Cin) \n) 



with dn = N — 1 at < and d„ = 1— ^i^,'^}^ at n > we see that the transformation 

" C(n+1) " C{n) — 



( I3.10p means 

G ^ C(0) exp (J2 ^J-k) G : exp {j^ d,r,i^^ 
k>i j 

for the group- like element. 

Equation f l3.8p encodes all the PDE's of the KP hierarchy. They are obtained by 
expanding the l.h.s. in the Taylor series in t' — t and equating the coefficients to zero. 
Technically it is convenient to set — )■ — t[ — )■ ti + ai and expand in — )■ 0: 

res^ [r(t - a - [z-^]) r(t + a + [z-^]) e-2?(a,^)] 
res, e«(^-^'')(r(t - a) r(t + a)) e^^^C^.^) 



res. 



z-^hj{da){T{t - a) r(t + a)) ^ z%{-2s.) 



.i>o 



l>0 



= J]/i,(-2a)/i,+i(4)r(t-a)r(t + a) = 0. 

i>o 

In the second line, the shift by [z^^^] is represented by action of the exponentiated differ- 
ential operator 

e(a.,.-^) = 5^^9„^.. 



The last equality can be written as 



h,(-2ci)hj+,{dx) eS'>i"'^^^r(t - X) r(t + X) 

j>0 



0. 
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Using the symbols Di for the "Hirota derivatives" defined by 

P{D)f{t) ■ g{t) := P{dx){f{t + X)g{t - X)) 



x=o 



where P{D) - is any polynomial in Z^j, we can rewrite it in the form 

^ /ij(-2a)/ij+i(£)) e^'>i"'^'r(t) • r(t) = (3.12) 

i>o 

The first non-trivial equation contained here is 

{D{ + 2,DI-ADiD^)t -T = Q (3.13) 

(the KP equation in the bilinear form) or 

rriiii - AtiTui + 3 {ruf + 3rr22 - 3 (r2)^ - 4rri3 + AtiT^ = 0, (3.14) 

where r, := dt-T. The second derivative of this expression gives the KP equation in its 
standard form 

3M22 = (4M3 - 12UU1 - Mlll)i , 

where u = df_^ log(r). 

Equation (13. 9p encodes the PDE's of the hierarchy which is sometimes called 1- 
modified KP hierarchy. The first non-trivial equation is 

{Df - D2)rn+i ■ r„ = 



or 



dt2 log = dl log(r„+irn) + (dt, log 

Try. V r„, / 



(3.15) 
(3.16) 



3.2.2 The bilinear identity for the 2DTL hierarchy 

For the 2DTL case the bilinear identity is written in the form 

^ (n + 1| e^+^^+VfcGe-^-^*-) \n) {n' - 1| e-^+^^+VfeCe-^-^*'-) \n') 

k 

= ^ + 1| e^+^^+^GV-fce-^-^*-) \n) {n' - 1| e-^+^''+'>G^le-^'''^'-^ \n') 



or 



res. 



res^ 



z-^ {n + 1| e'^+(*+V(^)G'e-^-(*-) \n) {n' - 1| e^+(*'+)^*(^)G'e--^-(*'-) \n') 



{n + 1| e^+^'+^Gij{z)e-^-^'-^ \n) {n' - l\e^+^''+^Gij*{z)e''^-^''-^ \n') 
Using formulas f l2.65p and f l2.66p . we transform this to 

^~i^at+-t'+,z) ^ ^(^)e-^+(t+)Ge--^-(*-) |n) (n' - l| ^*{z)e'^+'^^'+^Ge''^-^^'-^ \n') 



res. 



res. 



z e 



1 f(t_-t'_,z-l) 



(n+l| e-^+(*+)Ge--^-(*-V(2) \n) (n'-l| e-^+^^'+^Ge'^'^^'-^*{z) \n') 
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Next, bosonization formulas fl2.93p . fl2.97p give: 



res. 



res. 



. e^+(tV)Ge--^-(t-+W) |n'-l> 



We can write this as the bihnear identity 

f ^n-n'g€(t+-t;,.)^^(t+-[;.-i],t_)r„,(tV + [^-i],t'_)dz 

«/ Coo 



(3.17) 



^"-"V(*-*'-'^"'V„+i(t+,t_-[^])r„,_i(tV,t'_ + [^])dz 

Co 



vahd for any n, ra', t±, t'_j_. The same bihnear identity can be obtained for the r-functions 
constracted as matrix elements (13. 5p of group-like elements with non-zero charge. The 
contour Cq encircles all singularities of the function e^*^*"~*'-'^ -"z"""'. In particular, if 
only a finite number of times are non-zero, then it is a small contour around 0. Note 
that at t_ = t'_ and n > n' the r.h.s. vanishes and we get the bilinear identity for the 
MKP hierarchy for the function r„(t) = r„(t,t_), as it should be expected. In a similar 
way, setting t+ = t'_,_ and n < ri' — 2, we see that the l.h.s. vanishes while the r.h.s. gives, 
after the change z — )■ 1/z, the bilinear identity for the MKP hierarchy for the function 
r„(t_) = r„(t+, -t_). 

The first non-trivial equation contained in bilinear identity (I3.17P is 

i DiD^iTn ■ Tn + T^+lT^-l = (3.18) 

4- 1 '^ri 



or 



9,,9,_,logr„ = -^I^±i^. (3.19) 

T 

n 

Subtracting these equations at + 1 and n one gets the 2D Toda equation 

dtA-i^n = e'^""'^"-^ - e'^"+i-^" (3.20) 

for = log(r„+i/r„). 

Remark. There is a freedom to multiply the r-function of the 2DTL hierarchy by an 
exponent of any linear form of times with constant coefficients: 

Tn{t+,t^) ^Cexp(Con+ C'fe4)rn(t+, t_) (3.21) 
corresponding to the following transform of the group-like element: 
G ^ Cexp ^) Gexp (coQ + J2 

\k>0 J \ fc>0 

Clearly, this transformation preserves the form of the bilinear identity. 
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3.2.3 3-term bilinear equations 

Setting n' = n and t'^ = tk — ^(-^f^ + + ^3'') (13 .Tp . we see that the essential 
singularity at 00 splits into 3 simple poles at zi, Z2,zs: 

g5(t-t',^) I 



Taking the residues, we arrive at the 3-term relation 

{Z2 - z^Yn (t - \z^^\) r„ (t - [^2^1] - [^3^1]) + (231) + (312) = 0, (3.22) 

where the last two terms are obtained from the first one by the cyclic permutations of 
the indices. Another possible choice is n' = n and = tfc + — \{^\^ + + 
then 

1 - — 

i{t-t\z) ^ £0 

(1-^)(1-^)(1-^)' 

V zi / V Z2'^ Zz' 

This leads to a slightly different (but equivalent) equation 
{z^ - z{){z2 - ^3)r„ (t - [z^^] - [z~^]) r„ (t - [^2"'] - [^3"']) + (231) + (312) = 0. (3.23) 

Setting n' = n — 1, t'l^ = tk — jizi'^ + Z2''), we see that the essential singularity at 00 
splits into simple poles at zi,Z2: 

\ Zi'y Z2' 

Besides the residues at these points, there is also a contribution from the residue at 00, 
so we obtain the 3-term relation 

Z2Tn+l (t - [^2^^]) r„ (t - [Z^^]) - ZiTn+i (t - [z^^]) r„ (t - [^2^^]) 

(3.24) 

+ {Zl - Z2)Tn+l{t)Tn (t - [2;^^] - [^2^^]) = 0. 

It can be formally regarded as a particular case of fl3.23p in the limit ^3 — )■ oo, ~^ 
0. The limit 2:3 —> 00 is smooth while the other one requires to assign a meaning to 
lim r„(t ~[zq ]). The form of the evolution multiplier, gests to set formally 

zo^O 

limr„(t±[zo"i]) = r„^i(t) (3.25) 

zo^O 

which actually holds if the function r„(t — [-2"^]) for arbitrary n can be analytically 
continued from a neighborhood of 00 to a neighborhood of the point z = Q and is regular 
there. If it is singular, as in the case of solutions relevant to quantum spin chains, (I3.25P 
should be substituted by a more general prescription which, however, also allows one to 
regard fl3.24p clS cl particular case of fl3.23p . 

Setting n' = n, t'^ = tj, - \ a-^, t'_^ = t_k - jb'' {k > 1) in (I3l7l) and taking the 
residues, we get the 3-term bilinear equation for the 2DTL hierarchy: 

Tn{t+ - [a-1], t_)r„(t+, t_ - [b]) - Tn{t+, t_)r„ (t+ - [a"!] , t_ - [b]) 

^ (3.26) 

= -r„+i(t+,t_-[6])r„_i(t+-[a ^],t_). 
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3.3 Schur function expansions and Pliicker coordinates 

3.3.1 Schur function expansions 

One may expand the MKP r-function in the Schur polynomials: 

r„(t) = ^CAHsA(t). (3.27) 

A 

The sum is over all Young diagrams including the empty one (for which C0(n) = t„(0)). 
The coefficients cx{n) ("Pliicker coordinates") can be determined from the formula (13. ip 
by inserting the complete set of states in between the operators 6*^+ and G and using 

(ETHD: 

r„(t) = J2 (^1 e^""^*^ 1^' ^) ^1 G \n) = ^{-lf^^\sx{t) (A, n\ G \n) , 

A A 

SO 

c^(n) = (-l)^W (A.nlGIn) 

(3.28) 

= (-l)^W {n\ . . . V';+„,(^j^n-/3rf(;,)-l • • • V'n-/3i-lG' \n) , 

where 6(A) is defined in f l2.15p . In other words, c\{n) are expansion coefficients of the 
state G\n): 

G \n) = ^CA(n) \\,n) . 

A 

For the KP r-function the Schur function expansion is 

r{t) = Y,cxSx{t), (3.29) 

A 

where c\ = ca(0) = (—1)^*^'*'^ (A, 0| G |0). For the 2DTL r-function the general form of the 
expansion is the double sum 

r„(t+, t_) = ^ cx^{n)sx{t+)s^{-t_) (3.30) 
\,fi 

where Cx^{n) = (^—lyW+bM (A,n| G n). The Schur function expansions of r-function 
are also discussed in [231 EH EH] . 

3.3.2 Restricted Schur function expansions 

If r(t) = X]A'^ASA(t) is a KP r-function, then the expansion 

rW(t)= J2 CASA(t) (3.31) 

eix)<N 

restricted to Young diagrams with not more than non-zero lines is also a KP r-function 
for all N > 0. This follows from the fermionic representation 



rW(t) = (Ole^+^^^P+jv^lO) 
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where P^^y is the operator introduced in ( 12.8ip . 

More generally, if r„(t+,t_) = XIa '^v(^)'^A(t+)sA(— 1_) is a 2DTL r-function, then 

rr'*^)(t+,t_)= J2 E cUn)sx{t^)s,{-t^) (3.32) 

£{X)<N+n e{^i)<M+n 

is also a 2DTL r-function for all and M (in particular, one of them can be equal to 
+00, so that the corresponding sum is unrestricted). This follows from the fermionic 
representation 

rW*0(t) = (n| e-^+(*+)Pl^ G P^m \n) , 

where P^jy and P~^j^ are the operators introduced in fl2.8ip . (Note that the operators Pl„ 
allow one to apply a similar restriction on the number of columns of the diagrams in the 
sums, while other projectors introduced in section 12.81 make it possible to impose more 
complicated restrictions.) In particular, if r„(t) = XIa '^A(^)sA(t) is a MKP r-function, 
then 

rrHt)= Yl CA(n).A(t) (3.33) 

e{X)<N+n 

is also a MKP r-function for all A^. This follows from the fermionic representation 
3.3.3 Determinant formulas for c\{n) of the Giambelli type 



Proposition 3.1 The following determinant formulas for the coefficients Cx{n) of the 
T-function hold: 

Cin) = (C0(n))-'^W+1 det C(„,,;3,)(ri) (3.34) 

i,j=l,...,d(\) 

(assuming that c%{n) ^ 0/ 

These are the Giambelli-like formulas where C(q,-|/3 )(?t,) are the Pliicker coordinates cor- 
responding to the hook diagrams {(y.i\Pj) = (a^ + 1, 1^^). In the last expression, we have 
taken into account that {n\ G \n) = r„(0) = C0(n). The proof consists in applying Wick's 
theorem in the form fl2^ to (1^:281) : 

cx(n) = (-l)^W((n|G|n))-'^W+V det (n|C+„^^„_,^._iG|n) 

ij = l,...,d(A) 

(3.35) 

= (C0(n))-^W+1 det C(.,,;,,)H- 

jj=l,...,d(A) 

Formulas of the Giambelli type for the expansion coefficients, or Pliicker coordinates 
of the r-function were given in [23] . Using the Jacobi identity for determinants, it is easy 
to see that they are equivalent to the 3-term bilinear relations for the coefficients c\{n) 
given in [2] (the Pliicker relations): 
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Here 

(a/lPji) = (ai, . . . , ^i, . . . , I8k,...,l3d), 



Another way to obtain this relation is to apply the Wick theorem to 

(A, n\ 1pn+ar'4^n+as'4^*n-p,-l'4^*n-p,-lG \n) . 

Equation ( 13.35^ allows one to obtain a useful representation for the state G where 
G is any group-like element with zero charge. Using f l3.35p . we can write: 

G \n) = cx{n) |A, n) 

X 

- V(c0(n))^"^ V det C(„,|/3,)(nX_^^_i . . . C^^^i^rH^^ • • • V^r>kvi 1^) 



d>l ai>Q2> ■•>ad>0 

/3l>^(2>...>/3d>0 



C0 (n) [l + XI ^^^'^^l^^ 5Z n 1/3,) (^X_/3i-l • • • V'^d • • • ^r^i 1^) 



d>l ai,a2. - ><»d>0 j = l 

/3l,/32,...,/3d>0 



C0(n)exp(' V |n) 
V ^-^ Cain) ^ / 



a,/3>0 ""^ 



(n 




l3-lG 


n) 


{n 


G 


n) 



a,/S>0 

(3.37) 

(it is assumed that {n\ G \n) ^ 0). 

3.3.4 Determinant formulas for Cx{n) of the Jacobi-Trudi type 

There are also determinant formulas of another type which connect the coefficients cx{n) 
with different n. They were obtained in |26j . 



Proposition 3.2 The coefficients C\{n) of the T-function obey the relations 

c\{n) = \ TT c^{n-k)\ det CA,-i+j(n - j + 1), (3.38) 

where Csiji) := c^s~i\o){n) = {n — l\ tp^j^_g_iG \n) are the expansion coefficients for one-row 
diagrams and it is assumed that 00(72) 7^ 0. An equivalent set of relations is 



c\(n 



n C0(n + k) det c^''-'+^{n + j - 1), (3.39) 

J--"- / i,j=i,...,\i 



, k=l 
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where d^iji) := C(^o\a-i){n) = {—ly {n + l\ipn-aG \n) are expansion coefficients for one- 
column diagrams. 

These relations are sometimes called quantum Jacobi-Trudi formulas. Note that the 
transformation Cx{n) — )■ C{n)c\{n) with arbitrary C{n) preserves the determinant for- 
mulas f l3.38p . ( I3.39p . Clearly, this freedom corresponds to the possibility of multiplying 
the r-function by any C{n), see f l3.10p . 

Let us outline the proof of f l3.38p . The main step is to recall Proposition 12.11 and 
apply formula f l2.13p for the state (A,n| in (13. 28 p . Using this we get 

cx{n) = (n-/3i-l| V;+A^^^^_(;3i+i) • • • \n) 

(recall that /3i + 1 = i{X)). Applying the Wick theorem in the form (I2.48p . we arrive at 
the desired result. The proof of (I3.39P is similar. 

Remark. Let A = (s°) be the rectangular Young diagram of height a and length s, and 
let c"(n) be the corresponding coefficients cx{n), then application of the Jacobi identitity 
for determinants to (I3.38P yields the following 3-term bilinear relation: 

c:{n)c:{n + 1) - c:^,{n)cU{n + 1) = c:-\ny+\n + 1). (3.40) 
3.4 Examples of r-functions 

In this section we consider the following familiar examples of r-functions: 

• Characters (Schur functions) 

• (Quasi)polynomial r-functions 

• Multi-soliton r-functions 

• Partition functions of matrix models 

and give their fermionic realization. 

3.4.1 Characters and (quasi) polynomial r-functions 

The Schur polynomials themselves are r-functions of the KP hierarchy. This follows from 
equation (I2.77p . 

r(t) = (0| e^-^(*) |A, 0) = (0| e^+(*)C,,-i • • • ViV^., ■■■^a, |0) = {-if^'^sxit) (3.41) 

and the fact that the product of ip- and '?/'*-operators in the correlator satisfies the BBC 
(I2.4ip . Note that one can represent the state |A,0) as action of a normally ordered 
exponent to the vacuum [27|: 

d(A) 
i=l 
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We call such r-functions characters because sx{t) with = ^tig^ is the character of a 
group element g G GL{N) in the representation with the highest weight A. If = ^ sti 
for g G GL{N\M), then s\{t) are super-characters. There are different possibilities to 
embed these KP r-functions into the MKP and the Toda hierarchies. One of them is to 
use the projectors fl230D . for example r„(t) = (n| e'^+(*) |A,0) (0|n) = 6n,oi-lf^^'^sxit). 

In fact the Schur polynomials are very special cases of more general (quasi)polynomial 
r-functions. In the case of MKP hierarchy these r-functions can be obtained as either 
mean values of group-like elements with non-zero charge between different vacua or as 
some special limiting cases of fl3.2l) . 

Fix two finite sets of distinct points pi G C, qj G C, and construct (finite) linear 
combinations of the operators ip{z) and their 2;-derivatives at the points Pi and of the 
operators ip*{z) and their z-derivatives at the points qj: 

*i(pO:=E«-^>W ' ^*{qj):=J2^jmdTr{z) • (3-42) 

m>0 m>0 ^ 

The product of N \E''s and M $*'s is a group-like element with charge N — M. Therefore, 

r„(t) = {n\ e-^+(*)^i(pi) . . . ^;v(P7v)$I(gi) • • • $ M(gAf) \n-N+M) (3.43) 

is a r-function of the MKP hierarchy. It is not difficult to see that it is a polynomial in all 
the variables n, tk multiplied by an exponential function of a linear combination of these 
variables (a quasi-polynomial). These solutions are sometimes called rational because 
ratios of such r-functions are rational functions. One can see that the Schur function 
expansion of the r-function fl3.43p is restricted to Young diagrams that do not contain 
the rectangular N x M diagram. In particular, at M = the summation is restricted to 
diagrams A such that £(A) < A^. 

Let us show how to obtain r-functions of the form fl3.43p by taking a singular limit 
of (13.21) . We illustrate the point by a very simple example. Consider the operator 
• gW (p)^*ir)m^ (with \r\ < \p\), where ip'{p) = dpip{p). It is a group-like element with zero 
charge for any (3 but at /3 = /3o = ((0| ip' {p)'ilj* {r) |0))~^ it is not invertible and equals 

.^^M'(v)r{r).^ ^ 1 ^ /3oV^'(p)V;*(r) - (5, (0|V^'(p)V^*(r) |0) = M{p)r{r). 

Therefore, r„(t) = (?t,| e'^+''*V(p)'^*(^) l'^) is a r-function of the MKP hierarchy. Then 
one can consider the limit r — )■ 0. The limit is singular and requires some care. From 

r"-V*(r)l^) =^r''~^~^^*\n) = ^ r"-'-Vri^) = |n - 1) + 0(r) 

I l<n-l 

we see that in order to get a well-defined limit, one should multiply the r-function by 
r"~^ before tending r — > (this is just a transformation of the form (13.101) ). Then we 
obtain the r-function of the form (13.431) : 

r„(t) = (n| e-^+(*V'(p) \n - 1) = (n - 1 + ^ Hfc/jp^-^e^^^'^). 

fc>i 

This example can be further generalized. Take operators "^iipi) of the form ( 13.42^ . 
As in the previous example, we can construct the group- like elements •e^^*^(P^)'^*(''^) • with 
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zero charge and consider the r-function (n| e-^+(t) ^e'^^*^^^^)'^*^'"^^ ; . . . •e^i'fiipiWiri) . |^^_ 
Again, we choose /3j = ((0| '^j{pj)ip*{rj) |0))~^, then each operator becomes the product 
Pj'^ j{pj)ip*{rj) (non-invertible). Thus the above r-function reduces to a r-function of 
the form 

/3i . . . (n| e'^+(*)^i(pi) . . . ^^(p^)^*(r^) . . . r^n) \n) . 

In order to implement the hmit Tj — )■ 0, we redefine Tj — )■ erj with e — )■ 0, then it is easy 
to verify that 

^(„_i)7v_ I jv(jv-i)^*^^^^) . . . il)*{en) \n) = (ri . . . rjv)-"+' AAr(r,) \n - N) + 0{e) (3.44) 
where A/v(rj) = det rf~^ = TTfrj — r,) is the Vandermonde determinant. This for- 

i,j=l,...,N * Al^ 
i>j 

mula says that in order to get a well-defined limit, one should multiply the r-function by 
^{n-i)N-^ N{N-i) i^Q^QYQ tending e — )■ 0. The rj-dependent factors (ri . . . risf)~'^~^^Aiy{ri) 
in the r.h.s. of fl3.44p are irrelevant because they can be eliminated by a transformation 
of the form (l3.1Up . We thus obtain the r-function 

r„(t) = {n\ e-'+(*)^i(pi) . . . ^n{pn) \n - N) (3.45) 

of the form (^M)- 

Remark. In this case the prescription f l3.25p is not directly applicable. The behavior 
of the function r(t ± [zq'^]) as zq can be found by using the bosonization formulas 
f l2.93p and moving the fermion operators ipi^o) or ip*{zo) to the very right position. In 
this way we obtain a more general version of the prescription f l3.25p : 

lim [(-zo)^'^r„(t T K'D] = ^«±i(t) (3-46) 

which is equally enough to deduce the bilinear equation fl3.24p from f l3.23p as a limiting 
case. 

3.4.2 Multi-soliton r- functions 

We begin with the most general soliton-like MKP r-function r„(t) = {n\ e'^+^'^^G \n), 
where G is given by 

N 



G = ^ exp(^ Akr{q^)HPk)) ^ (3.47) 

i,k=l 

with a nondegenerate matrix A and qi,pk € C such that \qi\ > \pj\- We have: 

Tn{t) = l+^Aik{n\ e^+il)*{qi)i){pk) 1^^) + ^ ^ AikAi>k'{n\ e-^+ip*{qi)ij*{qi>)ij{pk')ij{pk) \n)- 

i,k i,i',k,k' 

Reorganizing the summation and using (I2.18p . we get: 

N d 

^"(t) = i+E E det A^,^ det (^!^)rr(^)>.^-)-?(M^™). 

^ l<r,s<d l<r,s<d\qj —Pk^^^\qk ' 

l<ki<k2<...<k^<N 

(3.48) 



44 



In this expression one can recognize determinant of the sum of two matrices: 



r„(t) = det (J + AQ"'*) , (3.49) 

NxN ^ ^ ^ 



where the matrix Q"'* is 



(Q".t)^^ ^ ( ^y"e?(t-pO-€{t,..). (3.50) 



Qk -Pi vgfc 

The standard A^-sohton solutions of the MKP hierarchy are obtained in the case when 
the matrix A is diagonah 

N 

G = ^ exp(^afcZ^*(gfc)V^(pfc)) ^ (3.51) 

k=l 

In this case 

r„(t) = det (6,, + {p./qkre^^''''^-^^'''^A ■ (3.52) 

NxN y qi^-p. J 

More explicitly, this function has the form 

Tn{t) = 1 + J2e'' + J2 + E c.,Qfc9fce^'+''^+^'= + . . . , (3.53) 

i i<j i<j<k 

where 

Qi-Pi^Qi' ' (Pi - gi) (gi - ) ' 

The r-function of the A^-soliton solutions to the 2DTL hierarchy has a similar structure: 
r„(t+,t_) = e^*>i'=***-Hl + ^2^"^ + ^^J'^"'^''' + •••)' ^2-^^) 

where now 

^V, = pA"^€(t+,p,)-g(t+,g.)+g(t_,pri)-g(t.,gri)_ 
Qi-Pi^Qi^ 

There is another fermionic realization of the A^-soliton solutions to the MKP hierarchy. 
Let us introduce the following fermionic operators: 

^iiPi, Qi) := ^{qi) + bi^{pi). (3.55) 

They are group-like elements with charge +1. Therefore, 

r„(t) = {n\e'+^'^^iipuQi)---'^NiPN,qN)\n-N) (3.56) 

is a r-function which can be calculated using the Wick theorem in the form fl2.50p : 

x„(t) = det (n-j + l|e-^+(*)^i(pi,g,)|n-j) 

= ..det ^ (n-j + l|e-^+(*)(^(g,) + 6,^(p,))^:_^>-j + l) 

(3.57) 

= det {n-j + l\ (e^^^'^'Vlg.) + fc^e^^^'^'^Vfe)) \n-j + l) 

i,j=l,...,N 

= det (e«(*'^')gr^+6ie«(*'^^)pr^). 

i,j=l,...,N V ^« 1-1/ 
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One can show that this determinant representation of the A^-sohton solution is equivalent 
to (^Mi- 
lt should be noted that the multi-soliton solutions are degenerate cases of much 
more general algebro-geometric (or finite-gap) solutions associated with smooth complex 
algebraic curves (Riemann surfaces) according to the Krichever's construction. Their 
r-functions are essentially Riemann's theta-functions. The fermionic representation of 
the algebro-geometric r-functions was addressed in [28] . 

3.4.3 Partition functions of matrix models 

The unitary matrix model. As the first example, let us calculate the r-function 
corresponding to the group-like element P"*" (see (12.531) ): 

r^(t+,t_) = (Ar|e^+(*+)p+e-^-(*-) \N) . 

Here we follow [29l [30]. First of all, it is not difficult to see that r^f = at < and 
To = 1. For > 1 we have: 

TN = {0\ ^N-i ■■■^0 e-'+P+e-^-^o . . . ^jv-i |0) . 

Below we use the short hand notation introduced in ( ]2.80p . Using this notation, we can 
write 

TN = (0| e''+rN-ii-J+) ■ ■ ■ roi-J+)P^P^M-J~) ■ ■ ■ ^N-ii-J-)e~'- |0) . 

Equations (12.681) imply that the operators ipn{—J-), i'n{~'J+) ^^e formula for 
contain only positive modes and, therefore, commute with P'''. Moving one P"*" to the 
right and another one to the left, and using the properties mentioned above, we obtain 

TN = {o\rN-ii-J+)---ro{-j+)M-J-)---^N-i{-j-)\o) 

= det (0|V^;_i(-J+)^fe_i(-J_)|0) 

l<j,k<N ■' 

by the Wick theorem. The expectation value under the determinant can be represented 
as a contour integral as follows: 

(0| ^;(- j+)^,(- j_) |o) = J2 hait^)hi-t_) (o| |0) 

a,b>0 

The whole determinant can then be written as an A^-fold contour integral: 

= ]^ f ■ ■ ■ f IK'-. - '-'X-/' - ") n«*^-'-«'-'""i^ ■ (3.58) 
■ j<k 1=1 ' 
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In accordance with subsection 13.3.21 the Schur function expansion of this r-function co- 
incides with a restricted version of the expansion of the trivial r-function flA7p : 



rAr(t+,t_)= SA(t+)sA(t_). 

1{X)<N 

When t^k = —ik (the bar means complex conjugation), the expression ^{t+,z) — 
^{t-,l/z) is purely real for z on the unit circle and tn coincides with the partition 
function of the unitary random matrix model written in terms of the eigenvalues. In this 
form, it can be treated also as the partition function of a canonical ensemble of 2D 
Coulomb particles confined on a circle. 



Partition functions of normal and Hermitian random matrices. Let us fix an 

arbitrary measure d^{z) in the complex plane and consider the following group-like ele- 
ment: 



Go = : exp (^j^Mz)iJi{i/z)dfi{z) - E^^-^l) : 



(3.59) 



Here ^+{z) = "^^ipuZ^-, ^'+(-2) = ^^V'n^"" are truncated Fourier series containing only 

n>0 n>0 

positive modes. Obviously, Go commutes with P+. Expending the exponent into a series, 
one can represent Gq in a more explicit form: 



m /(p 

m=0 



iIj+^Zi) . . .i)+{Zm)P ipl{l/Zm)...ijlil/zi)dlJ,i...dfirn, (3.60) 



where dfij = dfi{zj) and the first term in the sum is P . 
Let us consider the expectation value 

r7v(t+, t_) = (A^l e-'+^^+^GoP+e-^-^*-) | A^) (3.61) 

and apply to it a chain of transformations similar to the ones made in the simpler case 
Gq = 1 previously considered. Again, r^f = at A^ < and tq = L For A^ > 1 we have: 

TN = {0\rN-i- ■ ■roe'-^GoP+e-'-tPo. . .i^N-i\0) 

= (0| e^+V>^„i(-J+) . . . ^o*(-^+)P+G'oP+^o(- J_) . . . ijN^i{-J-)e-'- |0) 

= (0| ^^„i(- J+) . . . roi--U)GoM-J-) ■ ■ ■ ^N-ii-J-) |0) . 
Substituting the explicit form of Gq, we get: 

X p-r+i^/z^) . . . r4^/z,)M-J-) ■ ■ ■ ^N-i{-j-) io). 

The next step is to notice that only the term with m = N contributes to the sum and 
all other terms vanish. Indeed, at m > A^ the state 

r+i^/z^) . . . r4^/z,)M-J~) ■ ■ ■ ^N-i{-j-) io) 
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is in fact the null state because the number of annihilation operators exceeds the number 
of creation operators while at m < the operator 

is in fact the null operator because P~ multiplied by the uncompensated positive ?/'-modes 
from the right gives (see f l2.55p ). Therefore, the expression simplifies to 



X p-r+{i/ zn) . . . r+{i/zi)M-J-) ■ ■ ■ i^N-i{-j-) io). 

Since there are as many annihilation operators to the right of as creation ones, the 
state that they produce from the vacuum is proportional to the vacuum state itself, i.e., 

V^;(l/^jv) . . . V^;(l/^i)^o(-J-) . . .^7V-l(-J-) |0) = |0) Cat, 

where the constant Cat is 

Cn = (o|^;(i/^^)...7/';(iM)Vo(-J-)...V'7v-i(-J-)|o) 

= det (0|^;(l/^,)^fc-i(-^-)|0). 

1<J,K<A' 

Because 

(oi r+{i/z,)'4^k-i{-j-) io) = 5^ 4 K{-t-) (oi io) 

aJ>0 

a>0 

the constant Cat is explicitly given by 

N 

Cn = /^N{zdX{e-^^'-'''\ (3.62) 
1=1 

where we use the convenient short-hand notation for the Vandermonde determinant: 

^N{zi) = det (zi~^) = Y\{zi - Zj). 

i>j 

Now, it remains to calculate 

(oi rN-A-J-^) ■ ■ ■ roi-j+)M^i) ■ ■ ■ m^n) io) = ^ det (o| rj-ii-j+)M^k) |o) 

N 



= det zi~'e^^'+'''^ = ANiz^)Y\e• 

l<j,k<N " ^ ' LA. 

1=1 

which can be done in a completely similar manner. Collecting everything together, we 
obtain the result: 

r^(t+,t_) = j;tj I |A^(^OI'n^^^*^'''^'^^*~'''^^'"(^')- (3.63) 
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Assume that djj,(z) = e'^^^'^^d^z is a smooth measure in the plane and t_k = —i^, then 
the expression C,(t^, z) — .^(t_, z) is purely real and the integral (I3.63P has a physical in- 
terpretation as the partition function of a canonical ensemble of N identical 2D Coulomb 
particles in the plane in the background potential W{z,z) = —U{z,z) + 27leJ2k^kZ^- 



1 r ^ 



det / z'-'z^-'e^^^'^^d'z. 



(3.64) 



It is proportional to the partition function of the ensemble of normal random N x N 
matrices $: Z^r oc J e^^^(*'*^\ with Zi being their eigenvalues. 

If the measure dfi is concentrated on a curve F C C, then the 2D integrals J^{. . .)d^z 
are reduced to ID integrals . .)\dz\ along F. This means that the 2D Coulomb particles 
are confined to the curve F. For particular choices of F the integral fl3.63p yields the 
partition functions of random matrix models of certain types in terms of eigenvalues. For 
example, if F is the real line, one obtains the partition function of Hermitian random 
matrices: 

AT 

r^(t+,t_) = — / {/\^{x,)f\{e^^'^-'-^'^^'^d^^{x,). (3.65) 

It is the r-function of the ID Toda chain because it depends on the differences tk — t-fc 
only. If F is the unit circle, the integral f l3.63p becomes identical to fl3.58p which is the 
partition function of unitary random matrices. 

The Hermitean two-matrix model. An expression (13.590 for the element Gq admits 
a generalization with a double integral with an arbitrary contour (see [311 [32] for details). 
For example 

Go = : exp ( / / e^yij+{x)r+{lly)d^i{x)d^2{y) - E^^'^l ) ^ (^.66) 

gives, through the same formula (13.610 . the partition function of the Hermitean two- 
matrix model: 

N 

r^(t+,t_) = — / A^(x,)A^(|/,) TTe^''^'+«(*+'^')-^(*-^')rf/xi(xOrf/i2(i//). (3.67) 

The Harish-Chandra-Itzykson-Zuber (HCIZ) matrix model. For an element Go 
similar to (I3.66P but with different choice of the integration contour and the interaction 
term, 



Go = : exp (^ff ^'''""'''MvmiM^^^ ~Vo^'^*) *' ^'"''^ 
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one obtains the partition function for the model of two coupled unitary random matrices: 

(t+, t_) = Ajv(wi)^^(^i) n ' 



(3.69) 

After the change of variables 

1 ^ 1 ^ 

^'^^") = ^-kib) = -^J2bi, k>0 

m=l m=l 

(the Miwa transform), this integral becomes equal to the partition function of the HCIZ 
matrix model: 

r^(t+(a),t_(6)) = c^^rr— - / [DU]exp{ctTUAU^B). (3.70) 



This integral is known as the Harish-Chandra-Itzykson-Zuber (HCIZ) integral 
It is a X unitary matrix integral with the Haar measure normalized by the condition 
InxnI^^I = 1, v4 and B are two diagonal matrices: A = diag (ai, 02, . . . , a„), B = 
diag (61, 62, • ■ ■ , bn) (see |Zj for details). The integral f l3.70p is known to be a solution to 
the 2DTL [271 [35]. 



3.4.4 Matrix models represented by diagonal group-like elements 

The simplest example of matrix models with diagonal group-like element Go is given by 
the model of unitary matrices fl3.58p . it corresponds to trivial Gq = 1. The operator 
representation fl3.6ip of such models is equivalent to the one suggested by Orlov et al 
In fact all models of normal matrices with axially symmetric measures dn, i.e., 
such that dfi{z, z) = dfidz]"^) belong to this class. For an axially symmetric measure, the 
bilinear form in the fermion operators in (I3.59P becomes diagonal: 

■^^ m,n>0 •^'^ n>0 

where 

Jc Jo 
(we assume that the measure is smooth with U{z, z) = U{\z\'^)), so in this case 




(3.71) 



and the r-function (13.6 ID does have the form (A^| e'^+e^e ^~ \N) with X = y^^ Xj ; ipjip* ■ 
In this case expansion fl3.30p is diagonal, i.e., cx^{n) = unless A = [36] • More precisely. 
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in our case we deal with a singular limit of the vacuum expectation value with Xj = log gj 
for j > and Xj — > +00 for all j < 0. Indeed, writing 

= n (1 + (^""^ - -11(1 + (^"""^ - 1)^1^^-) ' 

i>o i<o 

we see that the first product is equal to 11(1 + {gj — l)ipjip*) = Gq while the limit of the 

second one is the singular operator P+. 

Below we give three examples. One important example is U{z,z) = c|zp, then 

g,^ = -KC-'^-^nl = 7rc-"-^r(n + 1) (3.72) 

and Xn = —{n + 1) logc + logr(n + 1) (the common constant logvr is irrelevant). Note 
that the analytic continuation of this formula to negative values of n with the help of the 
gamma-function automatically implies the required singular limit X„ = +00 at n < 0. 
The r-function fl3.63p for this case has the following expansion in Schur functions [361 123] : 

N 

r^(t+,t_) = vr^c-^(^+i)/2nr(^)- E c-l^l(iV),.,(t+)s,(-t_). (3.73) 

k=l e{X)<N 

where the factor (A^)a is (A^)a := JJ(A^ + I - i){N + 2 - i) . . . {N + \i - i) (see also f lATj) 

i=l 

and Fig. 2 in Appendix A). 

The second example is the HCIZ model, where the element Gq given by fl3.68p has 
the same structure with g^ = c^/n\ which, up to an inessential ra-independent factor, is 
inverse to the one considered in the previous example (13.721) . so that the expansion of 
f l3.69p in Schur functions is as follows [131 137] : 

r^(t^,t_) = c^(^-i)/^nr4- E 7^^A(M.A(-t-). (3.74) 

k=l ^ ' e{X)<N ^ 

Our third example is the model of normal matrices with 

1 2 
U{z,z) = — (log|2/r|) . 

Here /3,r are parameters (r plays the role of a scale in the z-plane). The corresponding 
matrix model was introduced in [SS] (see also [35]). In this case we find: 

fOO 2 /"+00 2 

Jo J-00 

so, up to an inessential common factor. 

Go = exp ( 5^ ((/3 + 2 log r) n + ^ n'^^r^rn ] (3-75) 

\n>0 / 

51 



(see ( I3.7ip ). The Schur function expansion of the r-function (13.611) with this Gq reads: 

^An{N-1){2N-1)(^^P/2^N(N-1) ^ ^/3C;,/2/^^2^/3(7V+i)V^' 

e{x)<N 

where 



r^(t+,t_)=e^^(^-l)(2Ar-l)^^^/3/2)7V(7V-l) ^ ^/,C./2 j^^2^/3(7V+i)y ',,(t^),,(_t_), 

e{x)<N 

(3.76) 



1=1 ^ / 



3.4.5 Matrix models, cut-and-join-like operators and bosonization 

Sometimes for description of matrix models and related r-functions instead of using the 
group-like elements (13.591) it is more convenient to use other group-like elements. The 
reason is that the truncated Fourier series is not very convenient for bosonization, see 
section 12.111 



For example, according to Morozov and Shakirov [lO], the partition function of the 
Gaussian branch of the Hermitian matrix model (I3.65P (here we consider only one semi- 
infinite family of times t = — t_): 



1 /■ ^2 

TNit) = — / JA^{x.,)rl[e~'^^^(''^^Uxi (3.77) 

can be represented in terms of action of the cut-and-join-like operator on the trivial 
r-function: 

?)v(t) = w-^'e-^^-^w^, (3.78) 

where the operator 

Ty_2 = V ( ahtatbTT, +{a + h + 2)4+6+2 T^riSr 1 (^-^Q) 

belongs to the W^^^ algebra and, in accordance with (12.102p . ^ := wdw ■ More precisely. 



dto 

N 



Tn 



k=l 



SO that T7v(t) = for negative N. 

The boson-fermion correspondence allows us to rewrite (13.781) as a fermionic vacuum 
expectation value 

?W(t) = (iV|e^+Wexp(Vf2)|iV) 
with the group element given by exponential of bilinear fermionic operator 

Vr^2 = ^ f :(50W)':^^^ = -res.(z-^:^*(z)9,V(^):) =5ZA;(A;-l)V^fc^*_2 (3.80) 
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where we used, in particular, expansion f l2.92p . The matrix model r-function fl3.77p can 
be obtained with the help of the projection operator: 



where Gq is a diagonal operator of the form (13.711) with coefficients up to inessential 
factor coinciding with (13 .72 p . namely gn = \/2nT(n + 1). 

Similarly, the group-like element for the r-function (I3.76P can be rewritten in a more 
symmetric way: instead of (I3.75P one can use G = expiW^) 



= E ((/^ + 2 log r) n + I n') l^^rn : 

nez 

(3.81) 

= -res, (^z-':r{^) + 2 log r) zd, + ^{zd,?^ ^(z) :) 

since, as it follows from (I2.54p . GP+ = GqP^. For Q = 21og(r) + | this group-like 
element gives the r-function of the two-component KP hierarchy, namely 



T(t+,t_) 

where 



e^+(*+) exp (^^ W,^^ Q^" e^-^"*") |0) (3.82) 



K = -^^f: l^'dz = -res, (z'^ :^*(z) (^{zd^f + zd, + V^(z) :) (3.83) 



and 



L^ = 2-f: (d<P{z)Y :zdz = -Tes,(^z-':r{z) (^zO^ + ^) (3.84) 

is known to be equal to the generating function for double Hurwitz numbers [M| HH |12] . 
The operator (I3.83P is the fermionic counterpart of the famous cut-and-join operator 

Wo=Y^ iabtatbT;2 + {a + b)ta+b- 

a,b=0 ^ 



dta+b dtadtt 



(more precisely, the cut-and-join operator does not include derivatives with respect to to, 
which do not affect the result) so that (I3.82p is equal to application of this operator to 
the simple tau-function (]A7P : 



r(t+,t_) =exp (^Wo\ exp l-J^kht-k] = ^e^^^^/^Q 

^ \ fc>0 / A 



SA(t+)sA(-t, 



(3.85) 

A generalization of the Hurwitz r-function is given by the vacuum expectation value 
with the group element 



G = exp I ^ 

\fc=0 



(k + iy. 
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where 




It describes the stationary Gromov-Witten theory on CP^ relative to two poles jl^. 
2DTL r-functions of this type, their representations in terms of free fermions and their 
Schur function expansions were considered already in [27]. 

Integrable properties of those vacuum expectation values with respect to variables Xk 
are considered in the next section. 

3.5 Another type of time evolution 

Here we consider another type of multi-time evolution of group-like elements and as- 
sociated r-functions. This class of r-functions has appeared naturally in such diverse 
physical and mathematical contexts as Gromov-Witten theory on CP^ [H], supersym- 
metric gauge theories [13 US], the melting crystal model of topological strings |171I1H] and 
hydrodynamic description of Fermi gas in one spatial dimension [IHl ED]. Their natural 
fermionic realization requires introducing evolution operators other than c^^^*-*. 

3.5.1 The Hamiltonian evolution operator 

Let us introduce the bilinear operators 

Hk = $^ri'^:^nC: = -res,(z-i:^*(z)(^a,)'=^(z):), A; > l, (3.86) 

which we call Hamiltonians in contrast to the currents fl2.ti2p 

The operators Hk commute for all A; > 1. Introducing an infinite set of evolution param- 
eters T = (Ti,T2, . . .), we unify the Hamiltonians in the generating operator 

H{T) = Y,TkHu. (3.88) 

k>l 

In the physical interpretation. Hi is momentum of the system of free non-relativistic 
fermions and H2 is the energy operator. Accordingly, for purely imaginary Tj's X = iTi 
and T = iT2 are physical space and time variables |19]. One may also consider 

nGZ 

which is the charge operator Q f l2.24p . 

With the help of the simple identity e^^"^" = 1 + {e^ — l)^/'„ ■?/'*, an easy calculation 
shows that 
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and, therefore, 

(3.89) 

(cf. f l2.65p ). We call e'^'-'^^ the Hamiltonian evolution operator. Its adjoint action is 
thus diagonal on the fermionic modes ?/'„, ■?/'*. The action on il){z) and '4'*{^) is more 
complicated. However, the action of the special operator e^^^^ is simple. It shifts the 
argument of the series il){z), il)*{z): 

(3.90) 

which is in agreement with interpretation of X = iTi as the physical space variable for 
the system on a circle. 

From fl3.90p it is clear that the basis states (I2.12p are eigenstates for the Hamiltonian 
evolution operator. Namely, we have: 

JJe''^"'^^-^!* \X,n) = c„e^^(") |A,n) , (3.91) 

where 

{g-6_i-6_2™...-fe„^ n < 
1, n = (3.92) 

gbo+bi+.-.+bu-i^ n > 

is the solution to the difference equation Cn+i/cn = e^" and 

d{\) 

B\{n) = ^(6„+a, - &n-ft-i) = X](^n+A,-i - K^j)- (3.93) 
i=i j>i 

Similar formulas hold for the left action to the dual basis vectors. 

Remark. The construction of the Hamiltonian evolution operator admits a "g- 
deformation" related to the quantum torus algebra and to the integrable structure of 
the melting crystal model |48j . 

3.5.2 The r-function in T-variables 

Given any group-like element G (with zero charge), let us consider the expectation values 

(n|Ge^We-^-(*)|n) or {n\e'+^'^e''^^^G \n) . (3.94) 

At fixed T, they are KP r-functions in the variables t (and MKP r-functions in t,n). 
We are interested in the cases when they are r-functions in the variables T. First of all, 
it follows from formulas fl3.9ip and their duals that at t = or G = 1 these expectation 
values are, at each fixed n, exponents of a linear form in the variables T and they thus 
represent the trivial r-function. 
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In what follows we set n = and consider the expectation value (0| e'^+^*^e^^'^-'G |0). 
Using (E3ZD and f lCT]) . (^M>, we have: 

G\0)=Y, cx |A, 0) = 5^ cj-'^ i<^.?</(-l/'=) 1^' 0) ' 

/3i>/32>...>/3d>0 

where A = (ai, . . . , ad\Pi, . . . , Pd), = {0\G |0) and 

e^m |A,0) = exp(5^T,5^(af - (-A-1)^)) |A,0) . 
fc>i 1=1 

Combining this with (I2.76P and using the Giambelli determinant formula for SA(t), we 
can compute the expectation value: 

d>0 ai>Q2> ■•>«d>0 
/3l>fe>...>/3^>0 



d 

X 



Now, following Orlov's idea [51], we are going to put t = a[tL'^-'^] = {aw^^, | | . . .) 
for any complex a, w and use equation (]A15P from Appendix A for the explicit value of 
the Schur function s f^^^j^^{a[w'^]) . We get: 



d>0 ai>ti2> ■•><id>0 
/3l>^(2>...>/3d>0 



i<ij<dVai + /9j + 1 



-|-r r(a + + l)r(l - g + A)e^''^^^^ 
^ A-j- ati;"*+ft+ir(ai + l)r(A + 2) r(a)r(-a) 



where r(x) is the gamma-function. Let us compare the r.h.s. with the multi-soliton 
r-function f l3.48p . In terms of the matrix 

_ C(,-i,fc,i) w-^-^'+ir(z + a)T{k - a) 
(OIGIO) uT{{)T{k + l)T{a)T{-a) 

the former can be written as 

^0|e-^+([--^l)e^(T)G|0) .-^ , f k 



V V det A^k. det (-. 

^-^ l<r,s<d l<r,s<dVi„ — H 



(OIGIO) ^ ^ l<r,s<d '"-^^ l<r,s<dVv-l + fc, 

l<ki<k2<...<k^ 



d 

X 

m=l j>l 



a 

Y\ exjp(y^^Tj{{im-iy - {-K 
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After the identification pi = i — 1, qi = —i, i > 1, this expression formally looks like the 
oo-soliton r-function fl3.48p . More precisely, we have: 



^0|e'^^+([--^l)e^(T)G|0) 




^exp(j] A,fcV^*(-A;)7A(^-l));; |0) 



(0|G|0) 



i,k>l 



The form of the r.h.s. explicitly proves that the l.h.s. is the r-function of the KP hierarchy 
in the T-variables. We do not know what is the most general choice of the times t such 
that the vacuum expectation value (0| e^+^^h"^^^G |0) is a r-function in T-variables for 
any group-like G. 

Remark. As is argued in [50], a more general form of the r-function in T-variables 

is 



This fact was established in [50] using some field-theoretical arguments. Note that at 
G' = 1 and T' = the statement reduces to the example above after the redefinition 
G — )■ le"''^^^^lG. However, we do not know how to prove the more general statement by 
the same method or by operator methods developed in [18] . 



where G, G' 




57 

































^2 = 5 
















= 3 
















/33 = 


= 1 






= 3 
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Figure 1: The Frobenius notation. The Young diagram A = (9,7,6,3,2,1,1) = 
(8,5,3|6,3,1). 

Appendix A: Young diagrams and Schur functions 

We use the notation of [9]. 

Partitions and diagrams. A partition A = (Ai, A2, . . • , A^) is a sequence of positive 
integer numbers Aj such that Ai > A2 > . . • > A^. The number I = i{X) is called the 
length of the partition. The partitions are naturally represented by Young diagrams. 
The Young diagram of A is a table whose j-th row (counting from the top) consists of Xj 

boxes. We will identify partitions with diagrams and will denote the diagram of A by the 

e 

same symbol A. The total number of boxes in the diagram A is |A| = Aj. The empty 

1=1 

diagram is denoted by 0. 

By A' we denote the transposed Young diagram which is obtained from A by reflection 
in the main diagonal. Namely, X'j is the height of the j-th column of A. 

Given a Young diagram A = (Ai, . . . , A^) with i = i{X) nonzero rows, let (a|/3) = 
(ai, . . . , ad\(3i, . . . , /3d) be the Frobenius notation for the diagram A. Here d = d{X) is 
the number of boxes in the main diagonal and = Xi — i, /3j = A^ — i. In other words, 
ai is the lenght of the part of the i-th row to the right from the main diagonal and Pi 
is the length of the part of the i-th column under the main diagonal (not counting the 
diagonal box). Clearly, ai > 02 > . . . > ttd > 0, /3i > /32 > . . . > /3d > 0. If A = (a|/3), 
then X' = 0\d). Note that 

d 

J2{a^ + Pi) + d= |A|. 

1=1 

A box a; G A has coordinates if it is in the i-th line (from the top) and j-th 
column (from the left). Let m G C. One can define the generalized Pochhammer symbol 
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u-2 










u—3 











Figure 2: The Young diagram A = (5, 4, 1, 1) = (4, 2|3, 0). The generahzed Pochhammer 
symbol {u)x is the product of the contents of all boxes. 



{u)x as 



(Al) 



{u)x^ H iu + j-t) = {u)x,{u-l)x,...{u-e{X) + l), 

where {u)k = u{u + 1) . . . {u + k — 1) is the usual Pochhammer symbol (see Fig. 1). In 
the Frobenius notation, we have: 



d{\) 



(A2) 



or, explicitly. 



For example. 



hk{t) 



fei+2fc2+...=A; 



hi{t) = h , 

/l2(t) = \tl + t2, 

hit) = |i? + tit2 + h , 



^4(t) = ^ti + Ul + l f\t2 + tlh + tA ^ 



(A3) 



Schur functions. Let t = {^1,^2,^3, . . .} be an infinite set of parameters. The Schur 
polynomials sxit) labeled by Young diagrams A can be defined by the determinant for- 
mula 

sx{t) = . det hx^_i+j{t), 
where the polynomials hj are defined with the help of the generating series 

exp 

(Z! ^kZ^) = 1 + ^l(t)-^ + ^2(t)^' + . . . 



24"! ' 2 2 I 2 

It is convenient to put /io(t) = 1, hk{t) = for /c < and S0(t) = 1. From the definition 
it follows that dt^hk{t) = /ifc_„(t). The functions hk are elementary Schur polynomials 
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in the sense that for 1-row diagrams A = (j) with j boxes S{j)(t) = hj{t). Equivalently, 
one can define 

5A(t) = . det eA^_i+,(t), (A4) 
where the polynomials ej are defined with the help of the generating series 

exp(^(-l)'=-4fc2^) = 1 + ei{t)z + e2{t)z^ + ... 

k>l 

It is clear that 

e,(t) = {-iyh,i-t). 

For 1-column diagrams A = (l-^) with j boxes S(p)(t) = ej{t). Equations ( 1A3I1 . (1A4I) are 
known as Jacobi-Trudi formulas. It can be proved [9j that the Schur polynomials form a 
basis in the space of symmetric functions of the variables Xi defined by ktk = Yli ^i- 
terms of the variables Xi {i = 1, . . . , N), 

deti<i.j<Ar(xf"^^'"^ ) 

deti<jj<Ar (a;. 
We note that 

.,(t) = (-l)l^lsv(-t). 

It immediately follows form the definition that if — > a'^t^, then /i^ — )■ a'^hk and s\ 
a^^^sx (the quasihomogeneity property). 

The following formulas for Schur functions for hook diagrams A = = (a + 1, l'^) 

directly follow from the Jacobi-Trudi formulas: 

a 
m=0 

(A5) 

= {-lf^hp-m{-t)K+m+l{t). 

m=0 

The Giambelli formulas represent the Schur function for an arbitrary diagram A = (c5|/3) 
as a determinant of the hook ones: 

SA(t) = det S(„^|;3,)(t). (A6) 

l<i,j<a(A) 

The Cauchy-Littlewood identity. We note the Cauchy-Littlewood identity 

^SA(t)sA(t') = exp(^H,4), (A7) 

A k>l 

where the sum is over all Young diagrams including the empty one. Writing it in the 
form 

sxiy)sxid) = exp VkOt^^ 

X k>l 
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Figure 3: The hook length h{2, 3) = 7. 



where d = {dt^^ ^dt^, ^dt^,, ■ ■ ■} and applying to s^(t), we get the relation 



sx{d)s^{t] 



(A8) 



t=o 



which reflects the or t honor mality of the Schur functions. 



Skew Schur functions. Let /i C A be two Young diagrams. The skew Schur functions 
(or the Schur functions for the skew diagram A \ /x) are defined as 



where the Littlewood-Richardson coefficients c^^ are determined by 



5/.(t)s.(t) = 5^C>A(t). 



One also has [9] 

and, as an easy consequence, 

There are generalized Jacobi-Trudi formulas for the skew Schur functions ^ 
■SAv(t) = det /iA,-^,-i+j(t) = det ex'-^r^i+jit). 



(A9) 



(AlO) 



Hook lengths. The hook length at x = is defined as 

h{x) = h{i,j) = Xi + X'j - i - j + 1 
(see Fig. 2). An important quantity is the product of all hook lengths of a diagram 

Hx = l[h{x). (All) 
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In particular, for a hook we have 

1 1 



Set t, = (1,0,0,...), [1] = then 



SA(t*) 



and, more generally. 



sx{u[w ^]) 



sMi]) 



Ma 



(A12) 



(A13) 



(A14) 



(see p] for the proofs). Using (lA14p . ( ]A12p and the Giambelli formula together with the 
Cauchy determinant formula 



det 



we get 



S(al/3)l«F 



i<i7<d ai + + l Uii=Mk + A + 1) ' 
^ ^(_i)ftr(w + a, + i)r(i-w + A)\ n •<,'(«, -«,')(/3.-/3.') 



^■^ V r(«. + i)r(A + i)r(M)r(i 



u 



«^|'ire/=iK+A + i) 



(A15) 

Following |l9], we show how to prove flA14p for hook diagrams (a;|/3) using the 
fermionic operators. We start with 

.(„l^)(t) = (-1)^+1 (Ole-'+W^Vi^-lO) 



(27ri) 



^ /o C 



(0|e^+^*V*(C)V^(^) |0), 



where the integrals are over small contours around 0. From this it is clear that 

(a + /3 + l).(.|/3)(t) = ^^^^ / - / ^z-\-^-\zd. + Cdc) (0| e'^^'^riCmz) |0) . 

[zni) Jo z Jo ^ 

The vacuum expectation value under the integral at t = i.e., tk = uw~^ /k is 



Note that 



- w-cY C 



w-cY C 



w — z J C — z 



-uwC{w - z) " {w - C) 



w — z J C ~ z 

factorizes, so the double integral becomes a product of two ordinary integrals: 

{a + (3 + l)s^a\i3)iu[w-'^]) 



{-Ifu r dz 



w 



2niz 



w 



-u-l 



2mC 



1 - 



c 

w 



;-l)/^(M)«+i(l-n)^ 



a\ (3\ w 



which is flA14p for the hook A = (a| 
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Appendix B: Proof of Proposition 12.6 

Let G be a solution of f l2.4ip consisting of parts with various charges: 

G = ^ Gk, 

k 

where the charge of Gk is k. Assume that for some p and q with p — q = a > both Gp 
and Gq are non-trivial. Then from (12.411) it follows that 

^kGg ® tplGp = J2 Gq^k ® Gp^l (Bl) 

fcez fcez 

Let us define two functions: 

Z„(t+,t_) = (n+p|e^+(*+)Gpe^-(*-) \n), 

(B2) 

Wn{t+,t-) = (n+g|e^+(*+)G,e^-(*-) \n). 

To prove that at least one of Gp and Gg is identically equal to zero it is enough to show 
that for all m and m' (see Proposition 12. 2^ : 

Wm{t+,t.) Zm'{t+,t^) =0. (B3) 

From (IBll) in full analogy with (I3.17P for all n and n' we have 

/ ^"-"'-«e«(*+-*V.-)vr„(t+-[^-i],t_)Z„,(t; + [^-i],t'_)dz 

J oo 

(B4) 

= jf^"-"'e«(*-*'-^"')iy„+i(t+,t_-[z])z„,_i(t;,t'_+[z])rf^. 

We denote b = n — n'. Let us consider different values of b assuming that n is arbitrary. 
For 6 = a — 1 we put t'_^ = t_|_ and t'_ = t_. Then only the left hand side of (IB4p survives, 
and we get: 

iy„(t+,t_)Z„_,+i(t+,t_) =0. (B5) 

If —1 < 6 < a — 1 we first differentiate (1B4I) with respect to ta-6-i and then put t'^^ = t+ 
and t'_ = t . Again, only the left hand side survives, so we get: 

Wn{t+, t_) Zn-b{t+, t_) = - 1< 6 < a - L (B6) 

For 6 = — 1 we put t'^ = t+, t'_ = t_ so that only the right hand side of (1B4I) survives: 

iy„(t+,t_)Z„_i(t+,t_) =0. (B7) 



For — l<6<a — Iwe first differentiate (1B4I) with respect to t-b-i and then put t'^ = 1+ 
and t'_ = t_. Again, only the right hand side survives, so we get: 

Wn{t+, t_) Zn-b-2{t+, t_) = 0, -1< 6 < a - 1. (B8) 
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Therefore, combining flB5|) . flB6|) . flB7l) and flB8l) . we get 

iy„(t+,t_)Z„_fc(t+,t_) =0, 0<A;<a. (B9) 
Now assume that for some fixed n and < 

W„(t+,t_)Z„_fc(t+,t_) 7^0. 
Then we take a derivative of ( 1B4I) with respect to ta-fc-i and put t_|_ = t'_,_ and t_ = t'_: 

iy„(t+,t_)Z„_fe(t + ,t_) = I z'' (^—Wn+l{t+,t--[z\)\ Zn-k~i{t+,t^ + [z\)dz 

Jo yOta-k-i J 

so that W„+i(t+,t_) Z„_fc_i(t+, t_) 7^ 0. Thus, if — 2 < < 0, we have got a contradic- 
tion with (IBSp . if /c < —2, then we come to a contradiction after a finite number of steps. 
For k > a the argument is similar. 
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